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SUMMARY 


•  A  narrow  electron  beam  was  used-  to  selectively  determine  the 
component  densities  in  rarefied  helium-argon  flows  by  spectrally  analyzing 
the  beam's  luminescence..  After  a  systematic  study  of  the  electron  beam 
fluorescence  in  helium  and  argon  had  been  made  in  a  series  of  preliminary 
experiments,  point-byr point  measurements  were  made  along  and  off  the  axes  of 
underexpanded  free  jets  and  through  normal  shock  waves  produced  by  shock 
holders  placed  in  the  jets.  Various  argon  mole  fractions  were  used,  and 
the  jet  Reynolds  numbers  ranged  from  100  to  10,000.  The  argon  enrichment  found 
along  the  jet  axis  agrees  well  with  a  theory  recently  proposed  by  F.  Si  Sher¬ 
man.  Off  axis  tests  show  an  argon  deficiency  near  the  jet  boundary.  The 
experimental  shock  profiles  agree  only  qualitatively  with  one-dimensional 
theories.  Radial  diffusion  destroyed  the  otherwise  approximate  one-dimen¬ 
sionality  of  the  shock  waves,  resulting  in  a  significant  argon  enrichment 
inside  the  partially  stagnated  shock  holder. 
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NOTATION 

speed  of  sound 

correction  coefficient  in  the  viscosity  lav  (defined  by  Eq.  8) 

orifice  diameter 

binary  diffusion  coefficient 

Sherman's  diffusion  parameter  (defined  by  Eq.  7) 

mole  fraction  of  the  heavy  species 

electron  beam  current 

anode  current  of  photomultiplier 

atomic  or  molecular  mass 

Mach  number 

number  density  of  atoms  or  molecules 
gas  pressure 

radial  distance  measured  from  the  axis  of  a  free  jet 

radius  of  a  hemispherical  surface  centered  at  the  orifice 

Reynolds  number  based  on  the  orifice  diameter  and  the  stagna¬ 
tion  chamber  properties  of  the  gas 

Reynolds  number  based  on  the  radius  of  shock  curvature  and  the 
free  stream  properties  of  the  gas 

distance  along  a  streamline 

unit  vector  pointing  in  the  streamline  direction 

continuous  surface  everywhere  normal  to  the  streamlines 

Schmidt  number  based  on  the  stagnation  chamber  properties  of 
the  gas  (Sc0  =  V-j90  D12  ) 

gas  temperature  (absolute)- 

center  of  mass  velocity  of  an  element  of  gas 

diffusion  velocity  of  the  heavy  species  (defined  by  Eq.  3) 

volumetric  pumping  speed  of  the  tunnel 

distance  along  the  axis  of  a  free  jet  measured  from  the  center 
of  the  sonic  orifice 


location  of  hypothetical  source  of  the  free  jet  expansion 

axial  distance  to  the  center  of  the  Mach  disc 

distance  normal  to  shock  wave  measured  from  center  of,  shock 

thermal  diffusion  factor 
ratio  of  specific  heats  of  a  gas 
shock  thickness 

shock  thickness  based  on  the  maximum  slope  of  the  density 
profile 

ratio  of  the  pressure  behind  the  Mach  disc  to  the  ambient  pressure 
wave  length  of  light 

"Molecular  mean  free  path"  ahead  of  the  shock  (defined  by  Eq.l4) 
viscosity  of  the  gas 
mass  density  of  the  gas 

molecular  number  flux  fraction  of  the  heavy  species  (defined  by 
Eq.  5) 

stagnation  chamber  conditions 
sonic  conditions 

test  section  conditions  outside  of  the.  free  jet 
argon 
helium 
the  free  jet 

refers  to  the  heavy  species  > 
refers  to  the  light  species  . 
the  shock  wave 

refers  to  the  free  stream  conditions  upstream  of  the  shock, 
refers  to  the  conditions  downstream  of  the  shock , 
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I.  INTRODUCTION 

In  the  general  field  of  experimental  supersonic  gas  dynamics 
there  occur  many  instances  In  which  gas  mixtures  are  subjected  to  severe 
accelerations  and  decelerations.  Under  these  conditions  the  resulting  gas 
dynamic  forces  are  able  to  support  large  pressure  gradients,  accompanied  by 
corresponding  temperature  and  density  gradients.  In  many  cases;  these  grad¬ 
ients  can  become  quite  significant  over  the  distance  of  a  molecular  mean  free 
path,  and  diffusion  effects,  in  particular  those  resulting  from  barodiffusion, 
can  no  longer  be  ignored  and  must  be  included  in  a  general  description  of  the 
flow  structure.  Typical  examples  are  gas  dynamic  shock  waves,  and  Prandtir 
Meyer  expansion  fans  in  a  rarefied  gas  mixture. 

Since  planetary  atmospheres,  in  which  aircraft  and  spacecraft 
have  to  operate,  consist  of  mixtures  of  various  gas  components,  and  since 
the  gases  expelled  by  propulsive  devices  are  generally  made  up  of  constituents 
of  widely  varying  molecular  masses*  the  diffusion  phenomena  are  also  of 
practical  interest  in  aerodynamic  problems.  For  example,  the  heat  transfer 
to  the  surface  of  a  re-entering  space  vehicle,  the  combustion  characteristics 
of  a  fuel-air  mixture  in  a  supersonic  flow,  or  the  thrust  characteristics  of 
a  reaction  jet  (.rocket  exhaust),  are  all  influenced  by  the  distribution  of 
the  gas  components  in  the  relevant  gas  mixture. 

The  present  report  concentrates  on  two  specific  flow  configura¬ 
tions;  the  gas  dynamic  shock  and  the  underexpanded  free  jet  issuing  from  a 
circular  sonic  orifice.  Binary  gas  mixtures  of  non-reacting  inert  compon¬ 
ents  were' used  throughout ;  Helium  and  argon  were  chosen,  because  these  gases 
are  readily  available,  and  because  their  atomic  masses  differ  by  a  factor  of 
ten.  As  will  be  shown  in  Section  II  a  large  mass  difference  of  the  atomic 
species  makes  the  barodiffusion  effects  correspondingly  large  and  easily 
detectable . 


Whereas  the  structure  of  the  underexpanded  free  jet  and  that  of 
the  one-dimensional  shock  wave  in  pure  gases  are  reasonably  well  understood, 
the  same  cannot  be  said  for  the  case  of  gas  mixtures.  Experimental  investi¬ 
gations  of  the  species  separation  in  the  free  jet  have  been  performed  in  the 
past  for  two  reasons.  One  is  the  prospect  of  separating  isotope  mixtures;  the 
other  the  possibility  of  generating  high  energy  molecular  beams.  In  the 
latter  case,  a  light  gas  is  seeded  with  a  small  amount  of  heavy  gas.  In  the 
process  of  expansion  the  heavy  molecules  reach  essentially  the  adiabatic 
escape  speed  of  the  light  gas  and,  thus,  gain  a  larger  directed  mass  velocity 
than  could  be  obtained  by  adiabatically  expanding  the  heavy  gas  alone. 
Furthermore,  the  heavy  molecules  remain  on  the  average  closer  to  the  jet 
axis  than  the  light  molecules.  By  introducing  a  suitably  designed  stripping 
baffle  into  the  flow,  the  jet  core,  which  is  enriched  in  the  high  energy 
molecules  of  the  heavy  type,  can  be  isolated  for  use  in  molecular  beam 
experiments.  The  arrangement,  consisting  of  a  skimmer  in  a  free  jet,  has 
also  been  proposed  for  separating  isotopes. 

In  order  to  determine  the  most  efficient  combination  of  jet 
parameters  and  stripping  baffle  design  and  location  for  the  purpose  of  iso¬ 
tope  separation  or  molecular  beam  generation,  previous  experimenters  used 
a  mass  spectrometer  to  determine  the  composition  of  the  gas  extracted  from 
various  regions  in  the  free  jet  by  means  of  a  sampling  probe  or  skimmer, 
with  the  hope  of  gaining  an  understanding  of  the  diffusive  separation  occurr¬ 
ing  in  the  free  jet  itself.  The  results  obtained  in  this  fashion  by  Becker, 
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Bier  and  Burghoff  (Kef.  l)  and  by  Waterman  and  Stern  (Refs.  2  and  3)>  who  com¬ 
pared  the  composition  in  the  jet  core  behind  the  skimmer  with  that  in  the  jet 
periphery,  gave  an  inconsistent  picture,  of  the  actual  free  jet  separation. 

This  is  so  partly  because  the  samples  were  taken  from  a  finite  cross-section 
of  the  , jet  over  which  the  gas  composition  was  not  constant.  The  measured 
composition,  therefore,  represented  an  average  value  in  this  case,  and  could 
not  easily  be  related  to  the  actual  composition  at  any  specific  point  in  the 
free  jet.  Chow  (Ref i  4),  who  used  a  small  diameter  sampling  probe,  found  the 
gas  enriched  in  the  heavy  species  at  all  points  in  the  free  jet.  This  contra¬ 
dicted  the  measurements  reported  in  Refs.  1  to  3*  It  is  now  generally  believed, 
and  the  present  work  lends  strong  support  to  this  view,  that  Chow's  separation 
results  were  probe  induced.  In  fact,  Chow  has  recently  shown  (Ref.  5)  that  a 
curved  shock*  standing  in  front  of  the  sampling  tube,  could  account  for  the 
results  obtained  in  his  earlier  experiments.  Reis  and  Penn  (Ref.  6)  came  to 
the  same  conclusion  after  performing  some  sampling  probe  tests  of  their  own. 

It  is  also  suspected  that  the  compositions  measured  behind  skimmers  in  Refs . 

1  to  3  were  influenced  by  the  skimmer  interacting  with  the  supersonic  flow. 

In  the  present  experiments  all  probe  induced  effects  were  avoided  by  analyz¬ 
ing  the  gaseous  luminescence,  induced  by  a  thin  17.5  kv,  lOOuamp  electron  beam, 
for  the  determination  of  the  component  densities  in  the  undisturbed  free  jet. 

Recently  there  have  appeared  in  the  literature  several  attempts 
to  give  a  quantitative  theoretical  account  of  the  species  separation  in  a  free 
jet.  Of  these  the  papers,  by  Zigan  (Ref.  7)  and  Sherman  (Ref.  8)  are  of  inter¬ 
est  here.  In  Zigan' s  work  the  results  are  given  in  integral  form  and  allow 
a  comparison  with  skimmer  experiments  (e.g.  Refs.l  to  3) •  The  agreement  is 
poor;  presumably  because  of  flow  disturbances  introduced  by  the  skimmer. 

The  results  of  the  present  experiments  can  most  readily  be  compared  with 
Sherman  (Ref.  8),  who  gives  an  expression  for  the  mole  fractions  along  the 
\  axis  of  symmetry  of  the  free  jet.  This  is  done  in  Section  VII  of  this  report. 

The  structure  of  the  one-dimensional  normal  shock  wave  in  a 
binary  gas  mixture  has  been  treated  theoretically  by  several  authors  (Refs. 

9  to  ll).  Both  continuum  flow  concepts  and  kinetic  theory  have  been  applied 
to  the  problem.  So  far  no  experimental  tests  have  been  reported  for  the 
shock  wave  structure  in  gas  mixtures. 

The  motivations  for  the  work  reported  here  may  be  summarized 

as  follows: 

(1)  A  non-ambiguous  measurement  of  the  diffusive  separation  in  under¬ 
expanded  free  jets  was  in  order  to  settle  the  argument  as  to  what 
species  separation  exists  in  the  undisturbed  free  jet. 

(2)  A  direct  test  of  Sherman's  theory  for  the  axial  variation  of  the 
component  mole  fractions  became  feasible  with  the  use  of  the  electron 
beam  probe. 

(3)  Direct  experimental  evidence  of  species  separation  inside  shock  waves 
could  be  obtained  and  could  be  compared  with  existing  theories. 

(4)  By  probing  shock  waves  produced  in  front  of  a  partially  choked  shock 
holder,  placed  in  the  supersonic  flow  of  the  free  jet,  some  additional 
light  could  be  shed  on  the  "probe  induced  separation  effects". 
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(5)  The  results  could  be  expected  to  have  some  practical  application  in 
the  design  of  gaseous  diffusion  facilities  for  isotope  separation 
and.  in  the  design  of  seeded  molecular  beam  devices. 

II.  THEORETICAL  IREDICTIONS 

For  the  isentropic  flow  of  a  diatomic  gas  (7  =  7/5)  the  Mach 
number  variation  along  the  axis  of  the  underexpanded  free  jet,  issuing  from  a 
sonic  orifice,  has  been  calculated  by  Owen  and  Thornhill  (Ref.  12),  using  the 
method  of  characteristics.  More  recently  similar  calculations  have  been  per¬ 
formed  for  specific  heat  ratios  of  5 /3>  7/5 ,  and  9/7  by  Wolff  and  have  been 
reported  by  Sherman  (Ref.  13).  These  solutions  are  universally  applicable  to 
free  jets  expanding  into  an  ambient  gas  at  non-zero  pressure  in  the  region 
bounded  by  the  orifice  and  the  first  wave  front,  which  registers  the  existence 
of  an  ambient  pressure.  This  wave  front  is  usually  well  defined  by  the  barrel 
shock,  formed  by  the  coalescing  Mach  waves  reflected  from  the  jet  boundary, 
and  by  the  Mach  disc.  The  Mach  disc  location  is  given  by  an  empirical  law 
(Ref  s  .14  and  15') . 

xm/d  =  (2/3)  (pypji  ,  u) 

where  x^  is  the  distance  from  the  orifice  to  the  Mach  disc,  D  is  the  orifice 
diameter,  and  pQ  and  pbo  are  the  stagnation  chamber  pressure  and  ambient  pressure, 
respectively.  A  sketch  of  the  general  flow  configuration  of  an  underexpanded 
free  jet  is  given  in  Fig.  2.1. 

Once  the  Mach  number  is  known  locally,  all  other  thermodynamic 
gas  properties  can  be  obtained  from  the  isentropic  relations  in  the  inviscid 
case.  These  show  that  for  distances  larger  than  about  two  orifice  diameters 
the  gas  density,  p  (or  the  molecular  number  density,  n),  falls  off  as 
l/(x-x0)2,  where  xQ  is  the  location  of  a  hypothetical  source  near  the  orifice, 
given  by  the  intersection  of  the  streamline  asymptotes  with  the  jet  axis  (see 
Fig.  2.1).  It  should  be  noted,  however,  that  this  variation  of  the  density  is 
not  dependent  on  the  flow  being  isentropic,  but  is  strictly  a  result  of  the  geo¬ 
metry  of  the  source-flow-like  expansion  and  the  fact,  that  at  x/D  =  2  the 
velocity  has  virtually  reached  its  maximum  value  of  the  escape  speed.  For  this 
reason  density  measurements  which  agree  with  the  isentropic  density  profile 
cannot  be  used  to  establish  whether  the  free  jet  flow  is  isentropic  or  not. 

In  his  treatment  of  the  underexpanded  free  jet  of  a  binary  gas 
mixture  Sherman  (Ref.  8)  combines  the  continuum  flow  conservation  equations 
(Navier-Stokes) ,  written  for  a  two  component  mixture,  with  the  binary  diffus¬ 
ion  equation  based  on  the  Chapman-Enskog  theory  of  diffusion.  The  latter 
equation  may  be  written  in  the  form 

m  m  r  /  m  -m  \  "1 

pl-l  =  pD12  Vinp-aV^nTj-  Vf  ,  (2) 

where  the  subscript  "1"  refers  to  the  heavier  species;  f  is  the  mole  fraction 
of  the  heavier  species;  m  =  fn^  +  (l-f)  m2  is  the  mean  molecular  mass;  D12 
and  CL  are  the  binary  diffusion  coefficient  and  the  thermal  diffusion  factor, 
respectively;  is  the  diffusion  velocity  of  the  heavy  species  defined  by 

▼i  =  -  2  >  (3) 
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where  u^  is  the  mass  velocity  of  the  heavy  species  and  u  is  the  center  of  mass 
velocity  of  an  element  of  the  gas  mixture. 

In  the  free  Jet  expansion  and  in  the  shock  transition  the  tem¬ 
perature  gradient  is  in  the  same  direction-  as  the  pressure  gradient.  Thermal 
diffusion  opposes  barodiffusion  under  these  circumstances.  An  order  of  magni¬ 
tude  analysis  of  Eq.  (2),  as  applied  to  the  free  Jet  or  the  shock  wave,  re¬ 
veals  that  barodiffusion  is  by  far  the  dominant  separating  agent.  In  baro¬ 
diffusion  the  heavy  molecules  diffuse  towards  the  high  pressure  regions.  In 
the  free  Jet,  because  of  the  streamline  curvature  near  the  orifice,  high  pressures 
occur  on  the  Jet  axis.  Radial  barodiffusion,  therefore,  causes  an  enrichment 
of  heavy  molecules  along  the  Jet  centerline.  Once  a  separation  between  the 
component  species  has  taken  place,  ordinary  diffusion  tends  to  homogenize  the 
mixture  again. 


In  order  to  obtain  a  solution  for  the  set  of  equations,  Sherman 
expands  all  variables  in  terms  of  inverse  powers  of  the  Reynolds  number,  based 
on  the  stagnation  conditions  and  the  orifice  diameter:  e.g. 

f  =  f^  +  f^/Reo  +  f^2VReo  +  •••  > 
where  the.  Jet  Reynolds  number*  is  defined  by  the  expression, 


1 

Re  =  p  a  D/|i  =  m  n  (7RT  )2  D/p  .  (4) 

o  o  o  /lo  o  o  '  o'  '  o 


In  the  zeroth  approximation  the  flow  is  isentropic,  and  the  mixture  remains 
homogeneous.  Neglecting  all  terras  of  second  and  higher  order  Sherman  then 
solves  for  the  first  order  perturbation  terms  in  the  mole  fraction,  f,  and 
the  molecular  number  flux  fraction, 

0  =  nx  u^.s/(n^u^  +  n^.s  ,  (5) 

along  the  axis  of  the  Jet  by  integrating  the  equations  along  the  central 
streamline  (_s  is  a  unit  vector  along  the  streamline) .  Using  the  expressions 
derived  in  this  manner  in  Ref.  8  the  change  in  the  mole  fraction  of  the  heavy 
component  along  the  axis  can  be  written: 
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In  this  relation  M  represents  the  Mach  number, 
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*  The  Reynolds  number  used  here  is  twice  that  used  by  Sherman. 


4 


4  ’ 


and 


C  =  T  /{i  T  «  1 

rv  r>  * 


(8) 


The  Schmidt  number  is  defined  as  ScQ  =  Pq/p0  >  and  the  starred  quantities 
are  those  at  sonic  conditions.  For  a  monatomic  gas  mixture  (7=  5/3)  Eq*  (6) 
reduces  to 
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The  difference  between  the  molecular  number  fraction  and  the 
number  flux  fraction  on  the  jet  axis  is  given  by 
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which  for  7=5/3  reduces  to 
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The  fact  that  f-0  does  not  vanish  is  a  consequence  of  the  gas  components 
having  different  mass  velocities  in  the  x-direction. 


Equations  (6),(9)5(lO),  and  (ll)  are  subject  to  the  following 
assumptions  and  limitations: 

(1)  Sherman's  analysis  has  been  based  on  the  assumption  of  large  Reynolds 
numbers,  so  that  second  order  terms  (of  order  Re0"2)  in  the  asymptotic  expan¬ 
sion  could  be  neglected.  This  is  equivalent  to  stating  that  the  free  jet  ex¬ 
pansion  of  the  gas  mixture  as  a  whole  is  reasonably  well  described  by  the 
isentropic  theory,  and  that  the  diffusion  effects,  as  well  as  the  viscosity 
and  heat  conduction  effects,  introduce  only  small  perturbations  around  the 
isentropic  solution.  Furthermore,  as  a  consequence  of  retaining  first  order 
perturbation  terms  only,  ordinary  concentration  diffusion  is  neglected. 

Because  the  local  Reynolds  number  per  unit  length  decreases  with  increasing 
x/D,  it  is  evident  that  this  analysis  cannot  be  valid  for  arbitrarily  large 
values  of  x/D. 

(2)  Since  relatively  few  experimental  .data  exist  for  the  temperature 
dependence  of  the  diffusion  coefficient,  the  product  pD^2  in  Eq.  (2)  has  been 
replaced  by  \i/ Sc.  In  the  subsequent  analysis  the  local  Schmidt  number  has  then 
been  considered  independent  of  temperature,  and  it  has  been  assumed  that  the 
temperature  dependence  of  the  diffusion  processes  are  adequately  taken  into 
account  by  the  temperature  variation  of  the  viscosity,  for  which  sufficient 
data  are  available.  In  this  connection,  a  linear  viscosity-temperature 
relation  (d/d0  =  CT/T0)  has  been  assumed  which  becomes  increasingly  inaccurate 
with  increasing  x/D  and  the  consequent  very  low  temperatures,  since  for  a  real 
gas  the  relation  is  better  described  by  a  power  law  dependence,  the  exponent 

of  which  is  less  than  unity. 


Reference  8  presents  an  explicit  numerical  solution  for  the  case 
of  diatomic  gases.  In  order  to  find  the  corresponding  result  for  a  monatomic 
gas  mixture, Eq.  (.9)  has  been  integrated  numerically  by  the  author.  For  the 
integration  Wolff's  calculated  Mach  number  curve  for  7=  5/3  was  used  from  one 
nozzle  diameter  out  to  infinity.  In  lieu  of  other  available  data  this  curve 
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was  extended  backwards  through  the  orifice  using  impact  probe  data  (see  Ref. 

15)  for  the  case  of  air.  This  apparently  inconsistent  patchwork,  can  be 
justified  on  two  accounts. 

(1)  The  isentropic  expansion  of  a  gas  up  to  about  Mach  2  is  primarily  a 
function  of  the  geometry  of  the  constraints  on  the  flow  and  is  only  a  weak 
function  of  7  . 

(2)  It  has  been  demonstrated  in  Ref.  8  that  the  calculated  downstream 
separation  is  almost  insensitive  to  the  shape  of  the  axial  Mach,  number  profile 
through  the  nozzle. 

The  calculated  results  are  shown  in  Fig..  2.2,  These  indicate 
that  for  a  monatomic  gas  mixture  the  species  separation  takes  place  more 
rapidly  and  reaches  a  slightly  higher  value  than  for  the  diatomic  case.  Thus, 
it  appears  that  for  a  helium-argon  mixture  almost  all  of  the  species  separa¬ 
tion  in  the  expansion  takes  place  over  an  initial  distance  of  about  three  ori¬ 
fice  diameters,  and  remains  relatively  constant  thereafter.  Furthermore,  the 
theory  predicts  an  inverse  proportionality  of  f  -  f0  to  the  flow  Reynolds  num¬ 
ber.  In  fact,  the  species  separation  will  be  insignificantly  small  for 
ReQ  >  5000  for  all  mixtures.  For  a  given  Reynolds  number,  on  the  other  hand, 
the  change  in  the  mole  fraction  of  argon  appears  to  be  a  maximum  for  a  mix¬ 
ture  consisting  of  about  equal  numbers  of  helium  and  argon  atoms.  This  is  so, 
because  the  parameter  E  varies  from  zero,  for  fQ  =  0,  to  a  maximum  value  of 
about  three  at  fQ  =  0.5>  ana  back  to  zero  again,  for  f0  =  1. 

In  order  to  compare  any  experimental  data  with  the  theoretically 
predicted  species  separation,  the  jet  Reynolds  number  and  the  parameters  E  and 
C  must  be  calculated  for  the  particular  jet  under  investigation.  In  this  connec¬ 
tion,  the  product  pQD  serves  as  a  convenient  parameter  for  defining  the  experi¬ 
mental  conditions,  since  for  any  given  gas  or  gas  mixture  this  number  is  pro¬ 
portional  to  the  jet  Reynolds  number.  To  facilitate  an  easy  conversion  of 
the  experimental  pQD  value  to  the  jet  Reynolds  number  and  to  Re0/EC,  the 
quantities  Re0/pQD  and  Re0/EC  pQD  have  been  calculated  as  a  function  of  fQ 
for  helium-argon  mixtures  at  a  stagnation  temperature  of  300°K  (see  Fig.  2.3 
for  these  curves).  In  deriving  these  curves  Wilke’s  formula  (Ref.  l6)  was 
used  to  estimate  the  viscosities  of  the  gas  mixtures  from  the  pure  gas  viscosi¬ 
ties.  The  latter  were  calculated  according  to  the  Chapman-Enskog  theory 
assuming  an  atomic  interaction  potential  of  the  Lennard-Jones  type  (see  e.g. 

Ref.  17). 


For  a  theoretical  treatment  of  the  shock  structure  in  binary 
gas  mixtures  the  reader  is  referred  to  the  papers  by  Sherman  (Ref.  9)* 

Fujimoto  (Ref.  10)  and  Oberai  (Ref.  ll) .  The  following  two  features,  common 
to  these  theories,  should  be  noted. 

(1)  Owing  to  the  fact  that  barodiffusion  tends  to  accelerate  the  heavy 
molecules  in  the  downstream  and  the  light  molecules  in  the  upstream  direction, 
the  light  species  compresses  ahead  of  the  heavier  species,  resulting  in  an 
enrichment  of  the  light  gas  throughout  most  of  the  shock. 

(2)  The  overall  thickness .of  the  shock  transition,  in  terms  of  the  up¬ 
stream  flow  properties,  is  larger  than  for  a  pure  gas  under  similar  conditions. 

One,  very  unexpected  result  is  peculiar  to  Sherman's  analysis, 
however,  and  is  a  consequence  of  the  application  of  the  Navier-Stokes  equations 
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to  the  shock  structure  problem.  For  a  helium- argon  mixture  containing  two 
per  cent  argon,  Sherman  finds  the  heavy  gas  accelerating  in  the  upstream 
portion  of  a  steady  flow  Mach  2  shock  wave  before  undergoing  the  expected 
deceleration.  Oberai,  who  uses  a  Moti-Smith  type  analysis s  also  calculates  the 
shock  profiles  for  the  case  of  very  small  mole  fractions  of  the  heavy  gas,  but 
does  not  anywhere  observe  this  peculiar  behaviour  of  the  heavy  species.  Since, 
by  necessity,  both  these  theoretical  methods  are  approximate  treatments  of  an 
inherently  complex  phenomenon,  it  is  not  possible  at  this  time  to  discard  one 
analysis'  in  favour  of  the  other  on  purely  theoretical  grounds.  The  answer 
will  have  to  be  provided  by  experiment.  In  fact,  the  initial  behaviour  of  the 
heavy  component  at  low  concentrations  may  turn  out  to  be  a  very  sensitive  test 
for  the  adequacy  of  any  theoretical  model. 


III.  THE  ELECTRON  BEAM  DENSITY  PROBE 


The  optical  electron  beam  density  probe  was  suggested  and  de¬ 
veloped  by  Schumacher  and  Gadamer  (Ref.  18)  and  was  first  applied  to  an 
actual  flow  investigation  by  Gadamer  (Ref.  19),  who  showed  that  it  can  be 
used  to  accurately  determine  the  local  densities  in  a  rarefied  gas  flow. 

The  technique  employs  a  narrow,  well-collimated,  high  energy  electron  beam, 
which  is  projected  at  right  angles  to  the  flow  through  the  region  of  inter¬ 
est.  Because  of  the  relatively  small  collision  cross-section  of  the  gas 
atoms  for  20  kv  electrons,  the  interaction  of  the  beam  with  the  gas  is  weak 
and  results  mainly  in  the  scattering  at  small  angles  of  some  of  the  primary 
electrons.  For  gas  densities  below  4  x  10!5  atoms  per  cnP  (l20p  Hg  at  room 
temperature)  the  beam  remains  well  defined  and  virtually  unattenuated  over 
c-  distances  of  the  order  of  10  cm.  At  the  electron  energy  used  many  of  the 

electron-atom  interactions  are  of  the  elastic  type.  There  occur,  however, 
sufficient  inelastic  collisions,  resulting  in  electronic  excitation  and  subse¬ 
quent  de-excitation,  to  give  rise  to  a  luminescence,  which  is  largely  confined 
to  the  beam  region  itself.  By  focusing  the  luminescence  onto  an  optical  slit, 
which  is  at  right  angles  to  the  beam,  the  light  issuing  from  a  small  elemental 
gas  volume  is  selected  and  can  subsequently  be  analyzed  by  photomultipliers. 
The  photocurrent  measured  is  a  direct  indication  of  the  molecular  number  den¬ 
sity  in  the  elemental  volume  under  consideration.  The  component  densities  in 
a  gas  mixture  can  be  measured  selectively  by  spectrally  separating  the  light 
emitted  by  the  gas  components. 

It  should  be  pointed  out,  that  at  gas  densities  corresponding 
to  a  gas  pressure  of  0.2  torr  at  room  temperature,  the  average  time  between 
atomic  collisions  (of  order  10~°  sec  for  most  gases)  becomes  comparable  to  the 
life  times  of  the  excited  states  of  the  atoms  (of  order  10"°  sec  for  allowed 
transitions).  Urder  these  conditions,  a  measurable  fraction  of  the  electron 
beam  excited  atoms*  collide  with  other  atoms  (ground  state  atoms),  before  they 
can  emit  the  excitation  energy  in  the  form  of  radiation.  In  this  connection, 
it  is  quite  unlikely,  that  in  the  collision  all  of  the  excitation  energy  (of 
order  10  ev)  is  converted  to  kinetic  energy  of  the  two  colliding  atoms.  On 
the  other  hand,  it  is  not  difficult  to  think  of  collision  processes,  in  which 
the  excited  atom  is  knocked  into  a  different  energy  level,  from  which  radia¬ 
tive  de-excitation  takes  place  in  a  spectral  region,  that  lies  outside  the 
spectral  range  observed  (e.g.  in  the  ultra-violet).  This  makes  the  observed 
intensity  of  the  electron  beam  luminescence  less  than  one  would  otherwise  pre¬ 
dict.  The  effect  is  generally  referred  to  as  "collisional  quenching"  of  the 
radiation  concerned. 
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Since  the  collision  frequency  and,  therefore,  the  collisional 
quenching  effects  are  temperature  dependent  as  well  as  density  dependent,  it 
is  almost  imperative  to  work  with  gas  densities,  at  which  collisional  quench¬ 
ing  can  be  neglected.  For  the  visible  spectrum  of  nitrogen  Gadamer  (Ref.  19) 
found  no  quenching  effects  up  to  0.1,  torr,  and  Muntz  and  Marsden  (Ref.  20)  found 
that  for  the  green  line  of  helium  (  A  5016  5)  the  quenching  effects  became 
important  only  above  0.25  torr, 


iv.  description  OF  APPARATUS 
4.1  The  Low -Density  Tunnel 

All  experiments  described  in  Sections  V,  VII,  and  VIII  of 
this  report  were  performed  in  the  UTIAS  steady  flow  low  density  tunnel,  the 
basic  design  of  which  has  been  described  by  Enkenhus  (Ref.  21).  The  pri¬ 
mary  pumping  system  of  the  tunnel  consists  of  nine  booster  type  oil  vapour 
diffusion  pumps,  which  have  a  total  maximum  pumping  capacity  of  about  10,000 
liters/sec  of  air  at  0.03  torr.  The  pumping  speed  is  somewhat  lower  for  argon, 
but  nearly  twice  that  value  for  helium..  Since  the  speed  of  sound  in  helium 
is  about  three  times  that  in  air,  however,  the  exit  diameter  of  a  helium 
nozzle  cannot  be  made  as  large  as  that  of  an  air  nozzle  (producing  the  same 
Mach  number  jet  at  the  same  jet  pressure)  in  order  to  not  exceed  the  tunnel's 
pumping  capacity.  In  the  present  tests  the  mass  flow  of  helium  was  also  limit¬ 
ed  by  the  mechanical  backing  pumps,  which  have  a  combined  pumping  speed  of  a 
nominal  485  cubic  feet/minute  at  a  pressure  of  about  1  torr,  regardless  of  the 
gas  used. 


The  pumping  speed  limitations  of  a  low  density  tunnel  define 
the  maximum  size  of  the  expansion  field  of  an  underexpanded  free  jet,  which 
can  be  produced  in  the  tunnel.  As  was  pointed  out  by  Ashkenas  and  Sherman 
(Ref.  15),  the  distance  from  the  sonic  orifice  to  the  Mach  disc  is # simply 
proportional  to  the  square  root  of  the  volumetric  pumping  speed,  V  ,  of  the 
tunnel, independent  of  the  orifice  diameter;  i.e. 


XM  = 


(7rt0)“^ 


7±l\  (7+D/M7-U  .1 
.  2  )  V 


(12) 


With  the  present  pumping  capacity  of  the  UTIAS  tunnel 
following  maximum  distances  from  the  orifice  to  the  Mach  disc  can  be 


the 

achieved: 


air,  Xj^  =  18  cm, 
argon,  x^  =  l8  cm, 
helium,  xM  =  13  cm. 


If  the  tunnel  pressure,  p^,  is  kept  constant  at  the  value  at 
which  the  volumetric  pumping  speed  is  a  maximum,  in  order  to  maximize  the 
size  of  the  test  region,  then  the  mass  flow  (which  is  proportional  to  p0D^) 
is  fixed.  In  this  case  the  jet  Reynolds  number  is  inversely  proportional  to  the 
orifice  diameter.  For  example,  if  the  orifice  diameter  is  decreased  from  Dq 
to  D2,  the  stagnation  chamber  pressure  has  to  be  raised  by  a  factor  of 
(Di/D2)^,  for  the  equilibrium  test  section  pressure,  p^,  to  remain  unchanged. 

The  new  jet  Reynolds  number  is  then  larger  by  a  factor  of  Di/l>2*  Furthermore, 
since  x^  is  unchanged,  x^/D  has  increased,  and  the  maximum  Mach  number  in  the 
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jet  (just  in  front,  of  the  Mach  disc)  is  correspondingly  larger.  Therefore, 
a  wide  variety  of  Mach  number  ranges  and  jet  Reynolds  numbers  can  easily  be 
achieved  by  using  a  suitable  selection  of  orifice  sizes. 

4.2  The  Electron  Gun 

The  electron  gun  assembly,  complete  with  diffusion  pumps^aynd 
gauges,  is  mounted  on  a  traversing  tables  which  makes  a  sliding  peal  with  the 
top  of  the  test  section  of  the  UTIAS  low  density  tunnel.  Ah  externdi  tra¬ 
versing  mechanism  allows  the  electron  beam,  which  passe?  through  the  test 
section  from  top  to  bottom,  to  be  moved  parallel  to  itself  within  a  15  cm  by 
15  cm  square,  A  schematic  diagram  of  the  low  density  tunnel, with  the  electron 
gun  in  place  is  given  in  Fig.  4.1.  An  external  view^of  the  gun  and  the  tra¬ 
versing  table,  with  the  power  supplies  and  the  instrumentation  in  the  back¬ 
ground,  is  shown  in  Fig.  4.2.  ** 

The  electron  beam  is  produced  by  a  commercial "television  tube 
gun  (Philips,  Part  No.  902-564)  housed  in  the  inner  of  two  concentric  tubes, 
which  pass  through  the  traversing  table  into  the  test  section  of  the  tunnel, 
where  the  normal  operating  pressure  i sab e tween  10"  and  10"  torr,  To  pre¬ 
vent  deterioration  of  the  electron  gun  cathode  and  to  eliminate  scattering 
of  the  electron  beam  before  it  enters  the  test  section,  the  gun  chamber 
(formed  by  the  inner  stainless  steel  tube)  must  be  maintained  at  a  vacuum  of 
about  10"°  torr  during  operation.  This  is  achieved  by  passing  the:  electron 
beam  through  an  intermediate  pressure  chamber  before  ejecting  it  into  the 
test  section.  The  intermediate  pressure  chamber  is  formed  by  the  outer  tube, 
which  projects  somewhat  farther  into  the  tunnel  than  the  gun  chamber.  The 
gun  chamber  as  well  as  the  intermediate  chamber  are  constantly  evacuated  by 
two  separate  80  liters/sec  oil  vapour  diffusion  pumps  (Edwards,  Model  203B) . 

A  flexible  backing  line  connects  the  diffusion  pumps  to  a  two  stage  rotary 
pump.  The  pressure  in  the  intermediate  chamber  is  less  than  5  x  10”  torr 
under  normal  operating  conditions. 

From  the  gun  chamber  into  the  intermediate  chamber  and  from 
there  into  the  test  section  the  electron  beam  has  to  ^ass  through  two  aper¬ 
tures  of  low  gas  conductance.  These  restrict  the  inflow  of  gas  from  the  test 
section  and  confine  the  electron  beam  geometrically.  The  upper  aperture  is 
2  mm  in  diameter  and  12.5  mm  long.  The  lower  aperture  is  1  mm  in  diameter 
and  15  mm  long  and  is  followed  by  a  copper  block  with  a  \  mm  wide  slit  to  re¬ 
duce  the  beam  width  as  seen  by  the  optical  apparatus  for  a  good  spatial  re¬ 
solution.  In  both  cases  the  apertures  consist  of  cylindrical  holes  drilled 
through  copper  inserts,  which  can  easily  be  replaced.  The  lower  one  is 
electrically  insulated  from  ground  by  means  of  ceramic  washers  so  that  the 
electron  beam  current  intercepted  by  it  can  be  measured  with  a  micro-ammeter . 

Careful  alignment  and  focusing  of  the  beam  is  necessary  to  allow 
it  to  pass  through  both  apertures.  A  system  of  eight  electrostatic  deflection 
plates  mounted  below  the  gun  electrodes  has  proved  to  be  an  invaluable  aid 
for  bringing  out  the  beam  and  for  keeping  the  beam  current  constant  during 
the  tests.  The  voltage  between  the  plates  of  any  pair  can  be  varied  between 
-400  and  +400  volts.  The  deflecting  plates  are  4  cm  wide,  4  cm  long,  and 
lcm  apart.  With  the  aid  of  the  control  circuit  shown  in  Fig.  4.3  it  is 
possible  to  shift  the  beam  parallel  to  itself  by  a  single  adjustment  for  each 
of  the  two  horizontal  dimensions.  Similarly,  the  beam  can  be tilted  about  its 
point  of  incidence  at  the  upper  aperture  by  adjusting  a  single  control  dial. 
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Any  defocusing  effect  due  to  the  transverse  electric  fields  is  minimized  by 
always  keeping  the  electric  potentials  of  the  plates  of  eachjpair  equal  in 
magnitude,  but  opposite  in  polarity,  with  respect  to  ground,  so  that  the 
region  traversed  by  the  beam  is  effectively  at  ground  potential. 

The  alignment  of  an  originally  misaligned  beam  is  effected  as 
follows.  First  the  beam  is  shifted  parallel  to  itself  until  the  sum  of  the 
currents  intercepted  by  the  bottom  aperture  and  the  collecting  cup  is  a  maxi¬ 
mum.  This  ensures  that  the  beam  passes  freely  through  the  top  aperture. 

Then  the  beam  is  tilted  around  the  top  aperture  until  the  Current  collected 
by  the  Faraday  cup  is  a  maximum.  If,  at  this  point,  it  is  felt  that  the  beam 
alignment  can  still  be  improved  further,  small  trimming  adjustments  can  be 
made  by  repeating  the  above  sequence  of  operations.  With  this  rapidly  con¬ 
verging  procedure  a  perfect  alignment  of  the  electron  beam  can  be  achieved 
with  relative  ease.  As  a  final  step  the  electron  gun  lens  voltages  are  ad¬ 
justed  such  as  to  minimize  the  current  intercepted  by  the  bottom  aperture 
for  an  optimum  focusing  of  the  beam. 

The  cathode  potential  is  supplied  by  an  NJE  high  voltage, 
regulated,  DC  power  supply  and  may  be  varied  between  zero  and  minus  30  kilo 
volts.  The  anode  is  kept  at  ground  potential.  The  grid  and  lens  voltages 
are  obtained  from  a  voltage  divider  circuit,  which  allows  ample  adjustment 
of  the  lens  voltages  (up  to  15$  of  cathode  potential) .  A  rectifying  and  filter¬ 
ing  network,  which  floats  at  cathode  potential,  supplies  the  filament  current 
(see  Fig.  4.3).  A  separate  regulated  400  volt  DC  power  supply  (Kepco  ABC  425  M) 
furnishes  the  deflection  plate  voltages. 

A  mechanical  gun  tilting  mechanism  is  also  provided  in  case 
a  coarse  alignment  of  the  electron  beam  becomes  necessary  (see  Fig.  4.4). 

This  mechanism  consists  of  two  mutually  perpendicular  pairs  of  loops  made 
from  spring  steel.  These  are  attached  to  the  glass  pumping  stem  of  the  elec¬ 
tron  gun  and  bear  against  the  inner  wall  of  the  gun  chamber.  When  one  of  the 
loops  is  compressed  or  elongated  by  turning  a  threaded  rod,  the  gun  tilts 
around  the  anode  clips,  which  rest  in  grooves  machined  in  the  anode  block  of 
the  gun  chamber. 

To  shield  the  electron  beam  from  the  earth's  magnetic  field 
and  from  any  stray  fields  in  the  laboratory,  the  outer  tube  surrounding  the 
gun  chamber  is  made  of  mild  steel.  The  fact  that  the  tunnel  itself  is  con¬ 
structed  from  mild  steel  plate  also  helps  in  this  respect. 

The  bottom  aperture  in  the  electron  gun  assembly  can  be  closed 
by  a  remote-controlled  door  in  order  to  isolate  the  electron  gun  chambers  be¬ 
tween  tests,  when  the  tunnel  may  be  at  atmospheric  pressure.  The  electron 
gun  elements  are  protected  against  back- streaming  oil  vapour  (CVC  Convalex-10) 
from  the  diffusion  pumps  by  water  cooled  baffles  and  liquid  nitrogen  traps. 
Protective  automatic  valves  have  been  built  into  the  electron  gun  pumping 
system  to  guard  it  against  the  detrimental  effects  of  electric  power  failures. 
The  pumping  connections  and  protective  control  circuits  are  shown  in  Fig.  4.5. 

At  the  instant  the  electricity  to  the  pumps  is  turned  off, 
purposely  or  accidentally,  a  magnetic  isolation  valve  between  the  mechanical 
pump  and  the  diffusion  pumps  isolates  the  system  from  the  rotary  pump,  and  the 
rotary  pump  is  automatically  vented.  In  addition, the  gun  chambers  are  iso¬ 
lated  from  the  diffusion  pumps  by  two  self-closing  2"  butterfly  valves.  To 
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close  the  valves  the  valve  handles  are  turned  through  90°  by  falling 
weights,  which  are  released  when  the  electric  current  ceases  to  flow  through 
the  solenoids  that  normally  support  the  weights. 

When  the  electric  power  is  turned  on  again  the  venting  valves 
close  and  the  mechanical  pump  starts  up.  A  slip  clutch  is  provided  to  protect 
the  motor  from  starting  overloads.  A  time  delay  relay  keeps  the  magnetic  iso¬ 
lation  valve  closed  for  120  seconds,  by  which  time  the  pressure  in  the  pipe 
between  the  mechanical  pump  and  the  isolation  valve  has  dropped,  to  10" 3  torr 
or  less.  After  these  two  minutes  the  isolation  valve  opens ,  and  the  mechani¬ 
cal  pump  resumes  pumping  on  the  backing  line  and  the  diffusion  pumps.  The 
diffusion  pumps  remain  off  .until  the  starting  relays  are  activated  again  by 
a  push  button  control.  Similarly,  the  2"  butterfly  valves  have  to  be  opened 
by  hand,  and  the  solenoids  holding  the  weights  have  to  be  switched  on.  The 
above  devices  have  been  designed  to  keep  the  activated  electron  gun  cathode 
under  a  high  vacuum  at  all  times  and  to  keep  the  inner  gun  chamber  from  be¬ 
coming  contaminated  with  pump  oil. 

A  hot  cathode  ionization  gauge  and  a>  Pirani  gauge  are  used 
to  monitor  the  pressures  in  the  gun  chamber  and  in  the  backing  line,  respec¬ 
tively. 


4.3  The  Optical  Analyzer 

The  function  of  the  optical  analyzer  is  to  measure  the  inten¬ 
sity  of  the  helium  and  argon  radiation  issuing  from  a  small  gas  volume  in  the 
electron  beam.  In  accordance  with  this  requirement  the  instrument  must  have 
a  high  sensitivity  to  analyze  a  basically  very  weak  light  source,  and  it  must 
spectrally  discriminate  between  the  light  emissions  from  the  helium  and  the 
argon.  Furthermore,  only  light  which  originates  in  a  well  defined  region  of 
the  electron  beam  should  be  allowed  to  enter  the  analyzer. 

Good  sensitivity  is  obtained  by  using  an  optical  system  of  low 
f-number  in  conjunction  with  high  gain  photomultiplier  tubes,  selected  for  a 
maximum  signal- to-noise  ratio  in  the  relevant  spectral  regions.  Two  achromatic 
f/2  lenses  are  used  in  an  approximately  symmetric  configuration  to  form  a  real 
image  of  the  electron  beam.  In  the  image  plane  a  slit,  normal  to  the  electron 
beam,  is  used  to  spatially  define  the  height  of  the  section  at  the  electron 
beam,  from  which  light  is  received  for  analysis.  To  spectrally  separate  the 
helium  radiation  from  the  argon  radiation,  the  light  beam,  emerging  from  the  slit, 
issplit  into  two  parts,  each  of  which  is  transmitted  through  an  optical  filter 
and  onto  the  cathode  of  a  photomultiplier. 

A  narrow  band  interference  filter  [Baird-Atomic,  Inc.,  Type 
B-ll)  is  used  to  select  the  helium  line  at  A  5016A.  The  half  width  of  the 
filter's  transmission  band  is  12A.  The  intensity  of  the  transmitted  light  is 
then  measured  with  an  EMI-9502S  tube.  The  argon  signal  is  obtained  from  an 
EMI-9558A  tube  used  in  conjunction  with  a  red  filter  (Kodak,  Wratten  Filter 
No.  88A),  which  transmits  light  only  at  wave  lengths  above  7200a.  The 
EMI-9558A  tube  is  insensitive  to  light  above  85OOA.  Several  intense  Arl  lines 
exist  in  the  spectral  region  between  7200A  and  85OOA,  but  practically  no 
helium  lines  are  located  there.  The  spectral  response  curves  of  the  photo- 
multiplier-filter  combinations  are  given  in  Fig.  4.6.  These  curves  were  ob¬ 
tained  by  multiplying  the  quantum  efficiencies  of  the  photocathodes  with  the 
transmission  factors  for  the  filters,  as  given  by  the  manufacturers. 


11 


The  photomultiplier  tubes  and  the  optics  are  contained  in  a 
light  tight  housing  (constructed  from  t*"  thick  black  "Arboron"  panels )',  which 
is  rigidly  supported  from  the  outer  electron  gun  tube  inside  the  tunnel. 

This  arrangement  has  the  advantage  that  the  optics  remain  aligned  with  the 
electron  beam  as  the  beam  is  moved  from  point  to  point  in  the  flow  field  under 
study.  Furthermore,  the  f-number  of  the  instrument  can  be  kept  at'  a  favour¬ 
ably  low  value  by  putting  the  analyzer  inside  the  tunnel  close  to  the  beam. 

On  the  other  hand,  the  typical  gas  densities  (between  10“^  and  KT1  torr) 
inside  the  tunnel  are  such  as  to  favour  electrical  breakdown  between  bare 
conductors,  and  suitable  precautions  have  to  be  taken  to  prevent  a  glow  dis¬ 
charge  between  the  high  voltage  connections  leading  to  the  photomultipliers. 

Originally  the  light  was  directed  from  the  optical  slit  to 
the  photomultipliers  via  a  lens,  a  splitter  plate,  and  a  prism.  As  will  be 
discussed  in  the  following  paragraph,  this  set-up  proved  unsatisfactory  and 
was  later  revised.  The  original  design  is  shown  in  Fig.  4,7*  In  this  version 
the  electron  beam  image  is  of  the  same  size  as  the  actual  electron  beam. 

The  light  emerging  from  the  2mm  high  slit,  S,  is  bundled  into  a  parallel 
beam  by  the  lens,  C.  A  dichroic  splitter  plate  (Bausch  &  Lomb,  Inc.,  Type 
45-I-600),  D,  at  an  angle  of  45°  with  respect  to  the  light  beam,  transmits 
light  of  wave  lengths  between  3500a  and  6000A  and  reflects  the  red  light 
above  6500A  toward  the  argon  detector  (red  filter  +  EMI-9558A).  A  right 
angled  prism,  P,  is  used  to  reflect  the  "green",  light  beam  toward  the 
helium  detector  (interference  filter  +  EMI-9502  s).  ~ 

Preliminary  tests  have  shown  that  this  optical  arrangement  is 
sensitive  to  small  changes  in  optical  alignment.  For  example,  a  small  lateral 
shift  of  the  electron  beam  of  |mm,  due  to  stray  magnetic  fields,  causes  the 
beam  of  light  to  shift  over  one  cathode  radius  (5mm)  at  the  EMI-9502  s  tube, 
with  a  resulting  attenuation  of  the  helium  signal  of  about  50  per  cent.  Ow¬ 
ing  to  the  large  photocathode  (44  mm  in  diameter)  of  the  EMI-9558A  tube  the 
argon  signal  is  affected  less  severely.  It  is  evident  that  this  arrangement 
is  not  very  satisfactory,  since  the  calibration  curves  for  such  a  system  are 
subject  to  change  according  to  the  space  and  time  variations  of  the  stray  mag¬ 
netic  fields  in  the  tunnel.  Of  the  results  discussed  in  this  report  the  only 
measurements  made  with  this  optical  arrangement  are  the  preliminary  tests  dis¬ 
cussed  in  Section  5*6. 

To  overcome  the  difficulties  outlined  in  the  foregoing  paragraph, 
the  optical  components,  S,C,D,  and  P,  have  been  replaced  by  a  fiber  optical 
light  guide  (American  Optical  Co.).  This  light  guide  consists  of  two  separate 
bundles  of  0.0?6  mm  diameter  coated  glass  fibers.  For  each  bundle  the  en¬ 
trance  section  is  rectangular  (l.2  mm  x  15  mm)  and  forms  the  horizontal  slit. 

The  exit  sections  are  circular  and  rest  directly  on  the  optical  filters  (see 
Fig.  4.8).  The  dimensions  of  the  light  guide  are  such  that  the  front  face 
of  the  fiber  optics  is  located  about  3  cm  ahead  of  the  original  slit  plane, 
requiring  a  shift  of  the  electron  beam  image  plane  by  that  distance.  This 
makes  it  necessary  to  use  the  achromatic  lenses  in  a  slightly  non- symmetrical 
arrangement  resulting  in  an  image  "magnification"  factor  of  O.69.  Hence, 
the  ,15  mm  width  of  the  fiber  bundles  in  the  image  plane  corresponds  to  a 
width  of  22  mm  at  the  actual  electron  beam. 

An  "ultra- stable"  Keithley,  Model  M-242, power  supply  provides 
the  cathode  and  dynoae  voltages  for  the  multiplier  tubes.  Purely  resistive 
voltage  dividing  chains  (about  200  k-ohms  between  dynodes)  are  used  for  obtain¬ 
ing  the  dynode  voltages.  An  additional  voltage  attenuating  circuit, built  from 
wire-wound  precision  resistors,  allows  for  the  possibility  of  applying  differ- 


12 


ent  cathode  potentials  to  the  two  tubes,  using  a  common  power  supply  (see 
Fig.  4.9  for  the  wiring  diagram  of  the  photomultiplier  circuits).  The  anode 
currents  are  measured  with  the  aid  of  a  Keithley,  Model  4l4,  pico-ammeter. 

The  photomultiplier  bases  and  the  voltage  dividers  are  immersed  in  diffusion 
pump  oil  to  prevent  a  glow  discharge  between  the  tube  pins  or  across  the  dynode 
resistors.  The  latter  are  contained  in  a  separate  oil  filled  box.  This  pre¬ 
vents  the  heat  dissipated  by  these  resistors  from  raising  the  temperature  of 
the  photomultipliers  which,  in  turn,  would  cause  an  increase  in  the  dark 
currents  of  the  tubes.  The  tubes  are  screened  against  external  magnetic 
fields  by  mu-metal  shields. 


4.4  Nozzles  and  Shock  Holder 


The  underexpanded  free  jets  investigated  were  produced  by  lett¬ 
ing  the  gas  expand  from  the  "high"  pressure  stagnation  chamber  into  the  "low" 
pressure  test  section  through  a  sharp  edged  orifice.  The  orifice  diameter 
consisted  of  a  5  mm,  15  mm,  or  47.5  mm  diameter  hole  machined. in  a  brass 
end  plate,  which  was  fitted  onto  an  existing  converging  nozzle  of  larger 
exit  diameter.  The  5  mm  and  15  mm  orifice  plates  were  attached  to  a  conver¬ 
ging  nozzle  of  47.5  mm  exit  diameter.  The  47.5  mm  orifice  plate  was  fitted 
to  a  subsonic  nozzle  of  l4o  mm  diameter.  The  wall  thickness  at  the  orifice 
was  equal  to  0.03  times  the  orifice  diameter  in  all  cases. 

The  normal  steady- flow  shock  waves  studied  were  produced  by  means 
of  a  shock  holder  placed  in  the  supersonic  jet  flow.  The  shock  holder  is  de¬ 
picted  in  Fig.  4.10.  It  consists  of  a  cylindrical  tube,  50  mm  in  diameter, 
with  its  axis  parallel  to  the  flow.  The  downstream  end  of  this  tube  can  be 
partially  closed  by  means  of  a  conical  plug,  in  order  to  choke  the  flow  through 
the  tube  just  enough  to  produce  a  normal  shock  at  the  entrance  lip  of  the  shock 
holder.  Four  narrow  slits  are  cut  into  the  front  lip  of  the  shock  holder, 
so  that  electron  beam  measurements  can  be  made  throughout  the  entire  shock 
wave.  Each  slit  is  3  mm  wide  and  12.5  mm  long.  The  slits  at  the  top  and  the 
bottom  allow  for  the  electron  beam  to  be  moved  into  the  shock  holder,  whereas 
the  slits  in  the  side  make  it  possible  for  the  optical  analyzer  to  "see"  the 
beam,  when  it  is  probing  the  downstream  portion  of  the  shock  transition  in¬ 
side  the  shock  holder.  Because  the  width  of  the  slits  is  small  with  respect 
to  the  typical  shock  thickness  encountered  and  with  respect  to  the  shock 
holder ' diameter ,  the  influence  of  the  slits  on  the  shock  geometry  is  believed 
to  be  negligible. 

A  photographic  view  of  the  spatial  arrangement  in  the  test 
section  of  the  shock  holder,  the  nozzle  (5  mm),  the  electron  gun,  and  the  optical 
analyzer  is  given  in  Fig.  4.11.  The  photograph  shows  the  bottom  end  of  the 
electron  gun  (with  the  remote  controlled  door  closed).  The  support  arm, 
clamped  around  the  gun  tube,  holds  the  optical  analyzer  (on  the  right)  and 
the  bracket  for  the  Faraday  cup  (vertical  rod  on  the  left) .  A  separate  view 
of  the  optical  analyzer  (original  version)  is  furnished  in  Fig.  4.12. 
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V..  PRELIMINARY  TESTS 


5.1  Static  Calibration  Using  Pure  Gases 

The  photomultiplier  anode  currents,  Ip^,  were  calibrated 
against  helium  arid  argon  densities  by  admitting  pure  gases  into  the  tunnel  and 
by  measuring  the  test  section  pressure  with  a  McLeod  gauge  and  the  tempera¬ 
ture  with  a  thermocouple,  All  calibration  tests  were  performed  with  a  17.5  k 
volt  electron  beam  energy  and  a  100  p  amp  beam  current.  In  these  tests  the 
gases  were  leaked  into  the  tunnel  at  a  steady  rate  (monitored  by  a  flow  rater) 
through  the  auxiliary  test  section  at  point  A  (see  Fig.  4.1),  while  one  or  more 
diffusion  pumps  kept  the  tunnel  pressure  at  a  steady  value.  The  calibration 
curves  for  argon  and  helium,  obtained  in  this  manner,  are  shown  in  Fig.  5.1. 

The  slight  non-linearity  of  these  curves  is  caused  by  the  presence  of  a 
fluorescent  halo  around  the  electron  beam,  which  at  very  low  densities  ex¬ 
tends  beyond  the  field  width  observed  by  the  optical  analyzer.  The  scatter 
of  the  calibration  points  is  due  to  the  uneven  pumping  of  the  diffusion 
pumps,  causing  small  fluctuations  in  the  tunnel  pressure  and  a  corresponding 
uncertainty  in  the  McLeod  gauge  readings. 

The  optical  filters  used  were  not  completely  effective  in 
separating  the  argon  lines  from  the  helium  lines.  There  existed  some  degree 
of  "cross-talk"  between  the  filters  producing  small  error  signals,  which  had 
to  be  calibrated  as  well.  This  was  easily  done,  for  example,  by  recording 
the  small  argon  signal  obtained  when  calibrating  for  the  helium  radiation  with 
pure  helium  in  the  tunnel.  Likewise,  the  helium  error  signal,  due  to  the  pre¬ 
sence  of  argon,  was  obtained  simultaneously  with  the  argon  calibration.  These 
erroneous  signals  were  one  to  two  orders  of  magnitude  smaller  than  the  "true" 
signals.  The  calibration  curves  for  these  error  signals  are  given  in  Fig. 5*2. 

It  is  believed  that  the  two  Aril  lines  at  A  5009$  and  A5017$ 
were  responsible  for  the  helium  error  signal  in  the  presence  of  argon.  On  the 
other  hand,  the  "argon  signal"  observed  in  the  static  helium  calibration  can¬ 
not  be  attributed  entirely  to  helium  lines  between  A  7200$  and  A  8500$.  The 
only  helium  line  listed  in  this  range  in  the" MIT  Wavelength  Tables"  (Ref.  22) 
is  the  Hel  line  at  A  7281$.  Being  very  close  to  the  cut-off  wavelength  of  the 
red  filter,  the  intensity  of  this  line  would  have  to  be  an  order  of  magnitude 
above  that  given  in  Ref.  22  to  fully  account  for  the  observed  signal.  It 
appears,  therefore,  that  the  "argon  error  signal  due  to  the  presence  of  helium" 
was  caused  to  a  large  extent  by  impurity  gases  present  in  the  tunnel  during 
the  static  calibrations.  With  only  one  diffusion  pump  evacuating  the  tunnel, 
background  gas  pressures  as  "high"  as  10"3  torr  were  observed  in  the  tunnel. 
This  background  gas  was  the  result  of  small  leaks  and  of  desorption  of  gases 
from  the  tunnel  walls.  In  this  connection  it  should  be  realized,  however,  that 
any  atmospheric  argon  leaking  into  the  tunnel  cannot  account  for  the  observed 
argon  error  signal  either,  since  the  amounts  of  argon,  which  could  possibly 
have,  been  present  in  the  test  section,  are  at  least  two  orders  of  magnitude 
below  those  indicated  by  the  argon  error  signal.  Although  it  is  now  believed, 
that  most  of  the  argon  signal  observed  during  the  static  helium  calibrations 
was  caused  by  Ngj  no  specific  attempt  has  been  made  to  identify  the  excitation 
mechanisms,  or  even  the  spectral  bands,  involved.* 


*  The  band  or  bands,  responsible  for  the  observed  light  between  A  7200$  and 
A  8500$,  probably  belong  to  the  first  positive  system  of  N2. 
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This  type  of  argon  error  signal  could  be  expected  to  be 
practically  absent  in  the  free  jet  flows,  because  the  level  of  impurity  gases 
was  significantly  lower  under  these  conditions  than  it  was  in  the  static 
tests.  In  the  free  jet  runs  9  diffusion  pumps  were  used  (as  compared  to  1  or 
2  in  the  static  calibrations)  to  pump  away  the  large  volume  flow  of  gas,  there¬ 
by  significantly  reducing  the  impurity  level  in  the  test  section.  In  addition, 
the  test  section  gas  could  not  penetrate  into  the  supersonic  jet.  Free  jet 
calibrations  of  the  error  signals  could  easily  be  made  by  observing  simultaneous¬ 
ly  the  helium  and  argon  signals ,  obtained  while  moving  the  electron  beam  along 
the  axes  of  pure  helium  and  pure  argon  jets.  In  a  pure  helium  jet,  for  example, 
the  local  helium  density  was  obtained  from  the  helium  signal  and  the  static 
calibration  curve  for  helium  (Fig.  The  argon  error  signal,  caused  by 

the  radiative  emission  of  the  helium  in  the  free  jet  at  the  indicated  helium 
density,  has  then  been  plotted  in  Fig.  5«2.  It  is  evident  from  Fig.  5*2  that 
the  argon  error  signal,  observed  in  the  pure  helium  jet,  is  very  much  smaller 
than  that  obtained  in  helium  containing  small  amounts  of  This  remaining 
small  argon  error  signal,  which  increases  linearly  with  helium  density,  can 
be  associated  with  the  Kel  line  at  A7281,3a  (3s  ^So.“%  lp^o).  For  the  pur¬ 
pose  of  reducing  the  free  jet  and  shock  wave  results,  presented  in  Sections 
VII  and  VIII  of  this  report,  the  free  jet  calibration  of  the  argon  error  sig¬ 
nal  has  been  used.  On  the  other  hand,  for  the  case  of  the  helium  error  sig¬ 
nal,  measurements  obtained  in  a  free  argon  jet  are  in  perfect  agreement  with 
the  static  calibrations  (Fig.  .5-2), 

When  the  electron  beam  was  near  or  inside  the  lip  of  the  shock 
holder  the  observed  beam  intensity  diminished,  because  the  shock  holder  sur¬ 
faces  intercepted  some  of  the  light  which  would  otherwise  have  entered  the 
lens  of  the  optical  analyzer.  Some  light  reflection  from  the  shock  holder  sur¬ 
faces  occurred  as  well,  increasing  the  amount  of  light  received  by  the  lens 
somewhat  above  that  reaching  the  lens  directly.  The  measurements  very  close 
to  the  nozzle  were  similarly  affected.  It  is  also  conceivable  that  these 
grounded  metallic  objects  could  have  actively  influenced  the  intensity  of  the 
fluorescent  halo  (e.g.  by  removing  metastable  and  ionized  atoms  from  the  beam 
surroundings,  or  by  absorbing  or  emitting  secondary  electrons,  or  by  absorbing 
or  reflecting  ultra-violet  photons,  which  are  involved  in  producing  the  helium 
halo).  The  changes  in  the  helium  and  argon  signals,  which  occurred  when  the 
electron  beam  was  moved  into  the  shock  holder,  were  calibrated  for  various  hel¬ 
ium  and  argon  densities  under  static  conditions.  Some  such  calibration  curves 
are  given  in  Fig.  5*3*  These  show  that  the  main  effect  was  a  signal  attenuation 
due  to  light  interception  by  the  shock  holder  lip.  The  solid  curve  in  Fig.  5«3 
gives  the  expected  signal  attenuation  for  an  ideal  point-source  of  light  locat¬ 
ed  on  the  axis  of  the  shock  holder  (assuming  non-reflecting  surfaces).  Because 
of  the  presence  cf  a  halo  around  the  electron  beam  the  experimental  attenuation 
factors  for  the  helium  and  argon  signals  are  density  dependent  and  lie  above 
the  theoretical  point-source  curve.**  Furthermore,  these  calibrations  had  to 
be  repeated  every  time  the  shock  holder  was  adjusted  or  moved,  because  the 
optical  alignment  of  the  slit  in  the  side  of  the  shock  holder  critically  in¬ 
fluenced  the  amount  of  light  received  by  the  optical  fibers  in  the  analyzer. 

The  signal  attenuation  for  electron  beam  positions  very  close  to  the 
nozzle  was  also  calibrated. _ _ 

*  It  is  shown  in  Section  5.4,  that  the  static  calibration  curves  can  be  used 
to  infer  gas  densities  in  a  high  speed  flow  without  introducing  a  significant 
error . 

**  A  calibration  of  this  effect  was  also  made  using  nitrogen  as  a  test  gas  and 
an  interference  filter  centered  around  the  Nj!)  band  at  A39T^A.  In  this  case  the 
signal  attenuation  followed  the  theoretical  point-source  curve  closely.  No 
halo  has  been  observed  for  this  Nj?  band. 


5.2  Electron  Beam  Excited  Spectra  in  Helium  and  Argon 

Before  these  calibration  curves,  obtained  with  pure  gases  under 
stagnant  conditions,  could  be  applied  with  any  confidence  to  measure  partial 
densities  in  gas  mixtures  flowing  at  supersonic  speeds,  the  following  ques¬ 
tions  had  to  be  dealt  with: 

(1)  What  are  the  life  times  of  the  excited  states  involved?  Are  these 
short  enough  (<10"“  sec),  so  that  the  excited  atoms  are  not  swept 
downstream,  and  are  not  lost  from  the  field  of  observation,  before 
they  radiate? 

(2)  How  do  density  gradients  affect  the  symmetry  of  the  fluorescent 
halo?  Are  accurate  density  measurements  still  possible  when  gradients 
are  present? 

(3)  Is  the  radiation  from  the  helium  influenced  by  the  presence  of  the 
argon  in  a  gas  mixture  (or  vice  versa),  giving  rise  to  erroneous 
measurements  of  the  species  concentrations? 


In  order  to  answer  these  questions,  it  was  necessary  to  have 
an  understanding  of  the  excitation  processes  involved.  In  particular  an  ade¬ 
quate  knowledge  of  the  nature  of  the  fluorescent  halo  was  required  in  this 
respect.  The  fluorescent  mechanisms  are,  therefore,  discussed  in  some  detail 
here.  In  the  following  discussions  the  word  "luminescence"  is  used  to  des¬ 
cribe  the  light  emission  in  the  electron  beam  region  resulting  from  the  pri¬ 
mary  excitation  of  the  gas  by  the  high  energy  electrons.  The  term 
"fluorescence"  is  reserved  for  the  halo  resulting  from  secondary  excitation  by 
photons  or  secondary  electrons. 


For  the  helium  line  at  A  5015. 7$  it  has  been  shown  conclusively 
by  the  work  of  Lees  and  Skinner  (Ref.  23)  and  Wolf  and  Maurer  (Ref.  24)  that 
the  fluorescent  halo  is  caused  almost  entirely  by  the  resonance  spreading  of  the 
ultra-violet  radiation  of  wave  length  537$  resulting  from  the  transition 
3p  Ip£  -  Is  -^Sq,  (Fig.  5-4).  For  every  thousand  3P  states  produced  by 
primary  excitation  in  the  electron  beam,  24  atoms  de-excite  to  the  metastable 
2s  !s0  state  with  the  emission  of  green  light  at  A  5016$,  whereas  976  atoms 
de-excite  down  to  the  ground  state  with  the  emission  of  ultra-violet  light 
at  A  537A.  The  mean  life  time  of  the  3p  state  is  only  1.74  x  10"9  sec 
(Ref.  25) .  When  a  neutral  helium  atom  in  the  vicinity  of  the  beam  absorbs 
a  A  537A  photon  this  atom,  is  raised  to  the  3p  state  and  may  again  emit 
green  light  at  A  5016$  with  a  0.024  probability.  Because  the  A  537$  photons 
have  a  relatively  short  mean  free  path,  the  radial  outward  flux  of  the  ultra¬ 
violet  light  is  rapidly  attenuated.  Consequently,  the  visual  intensity  of  the 
helium  fluorescence  falls  off  faster  than  with  the  inverse  distance  from  the 
beam,  as  would  be  appropriate  for  the  cylindrical  geometry  alone.  Before 
all  the  excitation  energy  of  the  original  3p  states  is  completely  converted 
to  A  5016$  radiation  and  metastable  2s  ^Sq  states,  the  average  ultra-violet 
photon  is  absorbed  and  re-emitted  40  times,  as  is  evident  from  the  ratio  of  the 
transition  probabilities  involved.  Because  the  re-emission  events  are  com¬ 
pletely  random  with  respect  to  direction,  however,  the  photon  does  not  diffuse 
more  than  a  few  optical  mean  free  paths  away  from  the  electron  beam  in  the 
process.  The.  fact,  that  the  helium  calibration  curve  (Fig.  5.1)  becomes 
linear  at  5  x  lO1^  atoms  per  cm^,  is  an  indication  that  at  higher  densities 
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essentially  all  the  A  5016$  photons  produced  reach  the  photomultiplier.  This 
means  that  the  ultra-violet  resonance  spreading  is  confined  to  a  radial  dis¬ 
tance  of  about  1  cm  at  an  atomic  number  density  of  5  x  10^  atoms  per  cm3. 
Consequently,  the  atomic  absorption  cross-section  for  this  ultra-violet  radia 
tion  is  of  the  order  of  lO"^  cm^.* 


With  respect  to  isolating  the  mechanisms  involved  in  the 
helium  fluorescence  the  work  of  Wolf  and  Maurer  (Ref.  24)  is  of  particular 
interest.  For  this  reason  a  short  outline  of  their  experiments  and  conclu¬ 
sions  is  presented.  The  vacuum  system  used  by  Wolf  and  Maurer  was  divided  in 
to  two  chambers  connected  by  a  small  aperture.  The  lower  chamber  was  tra¬ 
versed  by  a  high  intensity,  low  energy  electron  beam  (0.5  to  2.0  amp  and  35 
to  80  ev)  parallel  to  the  dividing  wall.  The  analysis  concentrated  on  the 
ray  of  fluorescence  penetrating  through  the  aperture  into  the  upper  chamber. 
The  helium  pressures  ranged  between  0.01  and  0.25  torr.  The  most  conclusive 
parts  of  this  study  are  itemized  hereunder. 

(1)  A  spectrographic  survey  revealed  that  most  of  the  observed  fluores¬ 
cence  consistedof  Hel  lines  due  to  transitions  from  the  np  ^P?  states  to  the 
2s  -*-Sq  state  (e.g.  A  3613$,  A  39&5$,  anc*  ^  5016$).  Other  weaker  lines,  how¬ 
ever,  were  observed  as  well  in  the  fluorescence,  e.g. 

A  4388$  (5d  1D2  -  2p  l-PJ) 

A  3889$  (3P  3lg}i}2‘  25  3S1} 

A  3820$ 

A  4026$ 

A  4472$ 

A  5876$ 

The  lines  in  this  second  group  became  more  intense  relative  to  the  lines  in 
the  first  group  with  increasing  helium  pressure. 

(2)  In  all  tests  a  strong  magnetic  field  of  1000  gauss  parallel  to  the 
electron  beam  forced  the  scattered  beam  electrons  and  the  secondary  electrons 
to  move  in  tight  spirals,  and  thereby  prevented  these  from  moving  away  from 
the  beam  and  into  the  upper  observation  chamber.  Furthermore,  an  electric 
field  of  330  volt/cm,  applied  normal  to  the  ray  of  fluorescence,  was  found 

to  have  no  influence  on  the  spatial  and  spectral  intensity  distribution  of 
the  fluorescence.  It  could,  thus,  be  concluded  that  secondary  electrons 
and  ions  played  at  most  a  very  minor  role  in  producing  the  helium  fluorescence 
in  this  configuration. 

(3)  A  quartz  plate,  which  absorbed  the  far  ultra-violet  but  not  the 
visible  light,  was  placed  across  the  observation  chamber.  In  this  case  no 
fluorescence  appeared  above  the  quartz  plate.  This  strongly  suggested  that 
all  the  observed  helium  fluorescence  owed  its  existence  to  the  ultra-violet 
resonance  lines  between  500$  and  600$. 


1-  (nd  3p  _  2p  3P°) 


(4)  The  intensities  of  all  lines  found  in  the  helium  fluorescence  in¬ 
creased  linearly  with  elecbron  beam  current,  indicating  that  a  fluorescence 
of  metastable  helium  atoms,  diffusing  away  from  the  electron  beam,  cou  not 

•ic  2 

*  Lees  and  Skinner  estimate  an  experimental  value  of  2  x  1Q~  '  cm  and 
derive  a  theoretical  cross-section  of  10"^3  cm^. 
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have  been  involved. 


The  observations,  described  above,  conclusively  showed  that 
the  strong  lines  (np  -  2s  -^Sq)  observed  in  the  helium  fluorescence  re¬ 
sulted  from  the  resonance  spreading  of  the  ultra-violet  transitions 
(np  -  Is  ^-Sq),  as  described  previously  for  the  line  at  A  50168.  The 
presence  in  the  fluorescence  of  the  other  "unexpected"  lines  (corresponding 
to  transitions  from  excited  levels  such  as  ^D,  3pOj  3d,  etc.)  could  only  be 
explained  by  assuming  that  the  responsible  terms  are  produced  in  the  "halo" 
by  collisions  of  the  second  kind  Of  excited  atoms  in  the  1P^- states  with  neu¬ 
tral  helium  atoms  in  the  ground  state,  e.g.^ 

He  Op'1^)  +  He  (is  hQ)  -» He  (is  *&Q)  +  He  (3d  3d). 

An  actual  spectrogram  of  the  helium  spectrum,  excited  by  an 
electron  beam,  is  presented  in  Fig.  5*5*  Here  the  electron  beam  energy  and 
current  are  17*5  kv  and  200  p  amp;  the  beam  width  is  \  mm.  The  helium 
pressure  is  about  0.1  torr,  and  the  gas  is  at  rest.  The  slit  of  the  f/l4 
spectrograph  (Hilger,  Model  E498)  is  normal  to  the  electron  beam  image. 

The  distance  between  the  two  vertical  lines  drawn  across  the  spectrogram 
corresponds  to  the  width  of  the  electron  beam.  In  the  spectrograms*  of 
Fig.  5.5  all  Hel  lines  between  A  38OOA  and  A  5900a  show  some  fluorescent 
spreading.  This  spreading  is  very  pronounced  for  the  lines  at  A  39&5a  and 
A  5016A.  A  relatively  weak  fluorescent  halo  is  observed  for  the  lines  at 


A  3820& 

(6d  3d 

-  2p  3po) 

A  38898 

(3p  3p° 

-  2s  3s  ) 

A  40268 

(5d  3d 

-  2p  3po) 

A  4l448 

(6d  XD 

-  2p  ip0) 

A  43888 

(5d  XD 

-  2p  ^P°) 

A  44728 

(4d  3d 

-  2p  3po) 

A  49228 

(4d  h 

-  2p  V) 

A  58768 

(3d  3d 

-  2p  3po) 

It  is, therefore, evident  that  for  a  17.5  kv  electron  beam  the 
characteristics  of  the  helium  fluorescence  are  very  similar,  if  not  identical, 
to  those  observed  by  Wolf  and  Maurer  for  an  80  ev  beam.  The  nitrogen  bands 
in  Fig.  5.5  are  caused  by  the  background  impurity  level  of  about  10"3  torr 
in  the  tunnel. 


The  electron  beam  excited  spectrum  of  argon  differs  markedly 
from  that  of  helium.  Whereas  for  the  case  of  helium  there^exist  practically 
no  Hell  lines  in  the  electron  beam  induced  luminescence,  a  substantial 
portion  of  the  radiation  from  the  electron  beam  excited  argon  results  from 
transitions  between  excited  levels  of  the  argon  ion.  Since  in  the  present 
experiments  only  Arl  lines  were  utilized,  the  following  discussion  will 
concentrate  on  the  spectrum  of  the  neutral  argon  atom. 


*  The  exposure  time  was  about  1  hour. 


18 


III* - — 


A  simplified  energy  level  diagram  for  neutral  argon  is  pre¬ 
sented  in  Fig.  5*6..  Only  those  levels  sure  shown  which  account  for  most  of 
the  observed  radiation  in  the  blue  and  red  parts  of  the  spectrum.  The  wave 
lengths  of  spectral  lines  resulting  from  transitions  from  the  3p^  '^P  and 
3p5  5p  configurations  to  the  levels  of  the  3P^  configuration  are  tabulated 
in  Fig.  5*7.  It  should  be  noted  that  for  neutral  argon  the  spectral  terms 
cannot  easily  be  designated  by  the  usual  symbols  established  for  atoms,  which 
exclusively  exhibit  LS-coupling  (Russell-Saunders)  between  the  orbit  and  spin 
angular  momenta  of  the  excited  electrons,  because  in.  many  instances  the 
electrons  of  the  argon  atom  exhibit  various  degrees  of  jj-coupling.  For  this 
reason  the  energy  levels  for  a  given  electronic  configuration  are  often  simply 
numbered  consecutively  in  the  order  of  increasing  excitation  energy.  The 
information  presented  in  Figs.  5.6  and  5 .7  has  been  derived  from  Refs.  22, 

26  and  27. 


A  spectrogram  showing  the  electron  beam  excited  argon  lines  in 
the  blue  and  violet  regions  of  the  spectrum  is  given  in  Fig.  5.8*.  This 
spectrogram  was  obtained  in  a  high  speed  argon  jet  with  a  1  mm  wide  electron 
beam  perpendicular  to  the  spectrograph  slit  and  the  flow.  It  shows  a  number 
of  Arl  and  Aril  lines,  many  of  them  exhibiting  a  halo-like  fluorescence  (e.g. 
the  Arl  line  at  A4200.7)>  which  extends  farther  on  the  downstream  side  than 
on  the  upstream  side.  Almost  all  the  Arl  lines  in  Fig.  5*8  are  the  result  of 
transitions  from  the  levels  of  the  electronic  3P^  5p  configuration  to  those 
of  the  3p5  4s  configuration.  These  lines  are  subject  to  the  same  excitation 
modes  and  transition  rules,  and  should,  therefore,  exhibit  the  same  charact¬ 
eristics,  as  the  red  lines  belonging  to  the  3P^  -» 3p^  4s  system  of  trans¬ 

itions.  The  characteristics  of  any  red  Arl  line  (contributing  to  the  argon 
signal  of  the  optical  analyzer),  can,  thus,  be  inferred  from  the  characteris¬ 
tics  of  a  corresponding  blue  line  in  Fig.  5.8.  Pairs  of  corresponding  lines 
can  easily  be  determined  from  the  table  in  Fig.  5*7.  Figure  5.8  is  presented 
here  primarily  because  of  this  close  correspondence,  and  because  of  the  fact, 
that  a  good  spectrogram  of  the  red  argon  lines,  resulting  from  electron  beam 
excitation,  was  not  obtained,  owing  to  a  lack  of  sensitivity  of  easily  available 
spectroscopic  films  in  this  region.  Nevertheless,  the  red  argon  lines  are 
significantly  more  intense,  on  an  absolute  scale,  than  the  blue  lines,  as  can 
easily  be  seen  from  Fig.  5.9A,  which  gives  a  spectrum  obtained  in  an  arc  heat¬ 
ed  argon  plasma  jet.  Figure  5*9®  shows  some  of  these  same  lines  obtained  from 
the  argon  luminescence  produced  by  an  electron  beam.  The  strongest  lines 
excited  by  the  electron  beam  in  this  spectral  region  appear  to  be  those  at 

A  7504$  , 

A  7515$  , 

A  8104$  , 

A  8115$  (strongest  line), 

A  84o8$  , 
and  A  8425$  . 


*  The  author  is  indebted  to  D.J.  Marsden,  formerly  at  UTIAS  and  presently 
at  the  University  of  Alberta,  for  this  spectrogram. 
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All  these  lines  are  the  result  of  transitions  from  the  elec¬ 
tronic  3p*  4pconfiguration  to  the  3p^  4s  configuration.  The  Einstein  transition 
probabilities  for  these  intense  red  argon  lines  are  all  about  5  x  10‘  sec"1 
(Ref.  28).  In  the  case  of  argon  the  even  terms  of  the  3p^  4p  configuration 
.cannot  be  excited  from  the  even  ground  state  by  resonance  absorption  of 
electric  dipole  radiation.  Outside  of  the  region  traversed  by  the  primary 
electrons  the  excited  3P^  4p  states  involved  can,  therefore,  only  be  pro¬ 
duced  by  secondary  electron  excitation  and/or  by  resonance  absorption  of  red 
light  by  excited  atoms  in  the  3p^  4s  configuration.  The  electronic  3P^  4s 
configuration  gives  rise  to  four  terms,  two  of  them  metastable  (see  Fig.  5.6). 

It  is  believed  that  light  absorption  and  re-emission  by  metastable  argon 
atoms  accounts  for  a  substantial  fraction  of  the  observed  intensity  in  the 
argon  halo  (some  evidence  for  this  is  presented  in  the  following  sections). 

Atoms  in  the  metastable  4s  1°  (3p°)  and  4s  3°  (3pg) states,  which  originate 
in  the  electron  beam,  diffuse  into  the  surrounding  ga?  where  they  can  absorb 
and  re-emit  red  light  of  wave  lengths  7635.11>  7723*76,  8014.79,8115.31a  and 

7724.21,  7948.18a,  respectively.  Since  this  type  of  fluorescence  depends  on 
the  diffusion  field  of  the  metastables,  the  red  argon  halo  should,  thus,  show 
a  definite  non- symmetry  in  a  high  speed  flow. 

It  is  also  conceivable  that  excited  argon  atoms  in  the  short 
lived  states,  4s  2°  (3p°)  and  4s  4°(lp°),  contribute  to  the  fluorescence, 
since  these  states  can  be  produced  outside  of  the  beam  from  neutral  argon 
atoms  in  the  ground  state  by  the  absorption  of  the  ultra-violet  resonance  lines 
of  wave  lengths  1066a  and  1048A.  During  the  short  life  times  of  1.05  x  10“tt 
sec  and  3.0  x  10"9  sec  (Ref.  29)  of  the  4s  2°  and  4s  4°  states  the  argon 
atoms  can  absorb  and  later  re-emit  the  red  argon  lines  of  wave  lengths 
7272.94,  7383.98,  7514.65,  8006.16,  8103.69,  8424. 65X  and  7503-87,  8264.52, 

8408.21,  8521.44a,  respectively.  The  part  of  the  fluorescent  halo,  which  is 
due  to  this  cause,  should  have  cylindrical  symmetry  in  a  uniform  argon  flow. 

The  fluorescent  mechanisms  discussed  above  will,  of  course,  also  apply  to 
the  corresponding  blue  argon  lines  (see  Fig.  5*7). 

Having  discussed  the  possible  fluorescent  mechanisms  for 
helium  and  argon,  it  is  now  feasible  to  predict  how  the  intensity  per  unit 
volume  of  the  fluorescence  should  vary  with  electron  beam  current  and  with 
gas  density.  For  this  purpose  let  1^  be  the  electron  beam  current,  n^  the 
number  density  of  atoms  of  the  relevant  gas  in  the  electron  beam,  and  n-f 
the  number  density  of  atoms  at  the  point  outside  of  the  beam  at  which  the 
fluorescent  intensity  is  to  be  considered.  In  cases  where  strong,, density 
gradients  exist  the  values  of  n^  and  np  may  well  be  different,  and  it  is 
important  under  these  conditions  (with  regard  to  the  spatial  resolution  of 
the  electron  beam  probe)  to  know.  ,  whether  the  fluorescence  observed  is 
proportional  to  nfc  or  np.  On  account  of  previous  arguments  the  gaseous 
luminescence  and  fluorescence  per  unit  volume  should  exhibit  the  following 
dependence  on  beam  current  and  gas  density. 

(1)  The  intensity  of  the  primary  beam  luminescence  per  unit  length  of 
beam  is  proportional  to  n^Ifo  for  both  gases. 

(2)  At  very  low  helium  densities,  when  the  ultra-violet  radiation  flux 
is  not  attenuated  appreciably  by  absorption,  the  intensity  of  the  green 
helium  fluorescence  per  unit  volume,  at  a  given  position  in  the  halo,  is 
approximately  proportional  to  n^n^I^. 


20 


(3)  For  the  argon,  line  at  A  8115$  (for  example)  the  local  fluorescence 
per  unit  volume,  due  to  absorption  and  re-emission  by  metastable  argon  atoms, 
is  approximately  proportional  to  n^I^,2, , since  the  radial  diffusion  flux  of 
metastables  and  the  radiation  flux  of  A  8115$  photons,  originating  in  the 
electron  beam,  are  both  separately  proportional  to  n^I^.  Similarly,  the  local 
fluorescence  due  to  ultra-violet  excited  argon  atoms  (e.g.  A  7504$)  can  be 
shown  to  be  proportional  to  n^n^I^2. 


(4)  If  it  is  assumed  that  the  secondary  electrons,  resulting  from  primary 
ionization  events  in  the  beam,  have  energies  in  excess  of  13  ev  for  argon  and 
in  excess  of  20  ev  for  helium,  then  the  fluorescence  per  unit  volume  due  to 
excitation  by  secondary  electrons  (at  a  given  point  within  the  range  of  the 
secondaries)  is  proportional  to  n^nfl^. 


(5)  If  the  secondary  electrons  have,  on  the  average,  energies  below  those 
cited,  above,  they  can  only  excite,  the  metastable  atoms,  and  the  resulting  local 
fluorescence  per  unit  volume  is  approximately  proportional  to  n^Ifo,  For  the 
case  of  argon  this  fluorescence  would  have  features,  which  are  very  similar  to 
those  discussed  for  the  first  mechanism  in  statement  no.  3« 


(6)  The  total  intensity  per  unit  length  of  the  beam  luminescence  plus  the 
fluorescence,  integrated  from  the  beam  center  to  infinity,  is  again  proportion¬ 
al  to  %I^,  since  it  is  a  direct  function  of  the  beam  energy  transferred  to  the 
gas  by  primary  electron-atom  collisions. 


To  obtain  a  photomultiplier  signal  which,  for  a  given  beam 
current,  increases  linearly  with  gas  density,  it  is  necessary  to  either  look 
at  the  beam  luminescence  alone  or,  otherwise,  include  all  of  the  fluorescent 
halo.  The  first  alternative  is  impossible  to  achieve  in  practice,  since 
there  is  always  some  fluorescence  along  the  line  of  sight.  Also,  the  signal 
would  be  a  function  of  the  electron  current  density  in  the  beam,  rather  than  of 
the  total  electron  beam  current.  The  second  alternative  can  be  approached 
closely  by  using  a  long  slit  at  right  angles  to  the  electron  beam,  as  was 
done  in  the  present  case.  In  this  configuration  the  spatial  resolution  in 
the  direction  normal  to  the  beam  is  still  equal  to  the  width  of  the  primary 
electron  beam,  because  the  total  light  output  is  proportional  to  nb  and  not  nf . 

As  a  third  alternative  a  spectral  line  could  be  chosen  which  is  known  not  to  occur 
in  the  halo,  but  which  is  present  in  the  primary  luminescence.  Since  none  of 
the  intense  helium  and  argon  lines  exhibits  this  property,  this  would  mean 
working  with  extremely  low  light  levels. 

p 

The ^dependence  of  the  local  fluorescence  on  n  (for  some  me¬ 
chanisms  even  on  n3)  causes  the  calibration  curves  in  Fig.  5*1  to  curve  upward 
at  the  low  density  end.  At  higher  density  levels  the  fluorescence  does  not 
reach  significantly  beyond  the  field  width  of  22  mm  observed  with  the  optical 
analyzer,  and  the  calibration  curves  approach  a  linear  density  dependence  in 
accordance  with  statement  no.  6,  above.  It  is.  interesting  to  note  that  accord¬ 
ing  to  statement  no.  3>  above,  the  local  argon  fluorescence  should  vary  as 
the  square  of  the  electron  beam  current.  A  somewhat  non-linear  behaviour  of 
the  argon  signal  with  changing  beam  current  was  indeed  observed.  Consequently, 
it  became  necessary  to  work  consistently  with  the  same  beam  current  for  both 
calibration  and  experiment.  In  contrast  with  this  behaviour,  the  helium  sig¬ 
nal  varied  linearly  with  beam  current  for  all  helium  densities,  indicating, 
in  agreement  with  Wolf  and  Maurer  (Ref.  24),  that  resonance  absorption  and  re¬ 
emission  by  atoms  in  excited  states  do  not  play  an  important  role  for  helium. 


5.3  Optical  Profiles  of  the  Electron  Beam  in  a  High  Speed  Gas  Flow 

In  order  to  study  the  characteristics  of  the  halo,  a  series  of 
tests  were  performed,  in  which  a  narrow  slit  (0.3  mm  wide),  parallel  to  the 
beam,  was  traversed  at  a  uniform  rate  across  the  electron  beam  image,  while 
the  photomultiplier  output  currents  were  simultaneously  displayed  on  an 
oscilloscope  screen.  In  ~this  fashion  the  profiles  of  the  electron  beam  indu¬ 
ced  light,  as  seen  through  the  optical  filters,  were  obtained  in  the  super¬ 
sonic  flow  of  helium  and  argon  jets.  Some  of  these  profiles  are  reproduced 
in  Fig.  5*10  showing  the  effect  of  pressure  gradients  and  flow  velocities 
normal  to  the  beam  on  the  symmetry  of  the  optical  radiation  field.  Figures 
5.10a  to  5*10e  give  the  intensity  profiles  of  the  argon  and  helium  radiation 
at  various  points  along  the  axes  of  underexpanded  helium  and  argon  jets.  The 
stagnation  chamber  pressure  and  temperature  are  25  torr  and  300 °K.  The  ori¬ 
fice  diameter  is  5  mm.  The  flow  is  from  right  to  left.  The  length  scales  have 
been  determined  from  the  ratio  of  the  oscilloscope  sweep  speed  to  the 
scanning  speed.  It  should  be  noted  that  under  normal  operation  the  field  of 
view  of  the  optical  analyzer  corresponds  to  2.5  times  that  covered  by  the 
oscilloscope  screen  in  Figs.  5*  10b  to  5*10e,  The  high  intensity  peaks, 
corresponding  to  about  a  1  mm  width,  must  be  ascribed  to  the  primary  lumine¬ 
scence  in  the  electron  beam.  The  lack  of  asymmetry  of  these  peaks  indicates 
a  mean  life  time  of  less  than  10~7  sec  for  the  upper  excited  states  involved, 
as  expected.  In  Fig.  5*10a  the  beam  luminescence  shows  a  considerable  spread¬ 
ing  due  to  small  angle  scattering  of  primary  electrons  by  argon  atoms.  At 
a  concentration  of  2  x  10^°  argon  atoms/cm3  (or  600  p  Hg  at  room  temperature) 
this  effect  can  be  expected  to  dominate  over  the  fluorescent  light  spreading, 
which  cannot  be  discerned  here.  At  an  equivalent  room  temperature  pressure 
of  145  p  Hg  (Fig.  5.10b),  however,  the  argon  halo  is  clearly  evident.  The 
higher  intensity  of  the  halo  on  the  downstream  side  is  caused  by  the  non-sym- 
metrical  diffusion  field  of  argon  metastables,  which  are  apparently  responsible 
for  most  of  the  fluorescence.  As  expected  from  the  predominantly  n2  depend¬ 
ence  of  the  argon  fluorescence  the  relative  local  intensity  in  the  argon  halo 
at  an  equivalent  room  temperature  pressure  of  34  P  Hg  (Fig.  5.10c)  is  only 
about  l/5  of  that  at  145  p  Hg.  Figures  5»10d  and  5*10e  indicate  that  the 
symmetry  of  the  green  helium  fluorescence  is  relatively  unaffected  by  the  gas 
velocities  and  density  gradients  in  the  free  jet.  It  is  also  evident  that  the 
ultra-violet  penetration  into  the  surrounding  gas  decreases  (i.e.  the  halo 
contracts)  with  increasing  helium  density. 

Figure  5»10f  shows  the  intensity  profile  observed  through  the 
helium  filter  in  an  argon  jet.  The.  light  is  believed  to  be  due  to  the  Aril 
lines  at  A  5009.35A  (4p  4po,2  -  4s  4*3/2)  and  A  5017.16a  (4p*  2F°,  -  3d  2D3/2) , 
which  are  close  enough  to  the  helium  line  at  A  5016a  to  pass  through  the  J 
helium  interference  filter.  The  "helium  error  signal"  due  to  the  presence  of 
argon  varies  linearly  with  argon  density  (see  Fig.  5.2).  There  exists,  there¬ 
fore,  no  halo  due  to  secondary  processes  for  the  above  two  Aril  lines.  The 
delayed  light  emission  evident  in  Fig.  5*10f  must  be  associated  with  a  rela¬ 
tively  long  life  time  of  one  of  the  two  upper  states  involved.  The  mean  life 
time  in  question  is  about  10"5  sec.  Since  the  4p  4P°/p  state  is  known  to  have 
a  mean  life  time  of  less  than  6  x  10"®  sec  (Ref.  28),'  the  long  life  time 
probably  belongs  to  the  4p'  2F°/  s^a^e* 


5.4  Axial  Density  Measurements  in  Pure  Gas  Jets 


The  calibrated  electron  beam  probe  was  used  to  determine  the 
density  along  the  axes  of  underexpanded  free  jets  of  pure  helium  and  pure 
argon.  Within  the  range  of  the  static  calibration  curves  the  measured  densi- 
ties  followed  the  theoretical  density  curve  quite  well  (Fig.  5.1l).  This  meant 
that  the  possible  error,  caused  by  the  halo  distortion  due  to  a  flowing  gas 
with  density  gradients,  could  be  neglected.  These  free  jet  runs  were  then 
used, to  extend  the  calibration  curves  to  higher  densities  (Fig.  5*12).  As 
is  evident  from  Fig.  5*11>  the  density  rise  due  to  the  Mach  disc  sets  in 
relatively  early,  even  though  the  Mach  disc  location  is  approximately  that  given 
by  the  accepted  empirical  law  (Eq.  l).  This  suggests  that  the  barrel  shock 
and  Mach  disc  are  quite  diffuse  for  the  Reynolds  number  range  covered.  A 
flow  visualization  study  of  the  free  jets  (Fig.  5*13)  substantiates  this.  The 
photographs  were  obtained  by  uniformly  traversing  the  electron  beam  in  a  plane 
defined  by  the  flow  axis  and  the  beam  itself  during  the  exposure  of  the  film. 

The  technique  has  been  described  by  the  author  in  Ref.  30. 

5.5  Shock  Profiles  in  Pure  Gases 


Since  normal  shock  waves  in  pure  gases  hav  *  in  the  past  been 
extensively  studied,  both  theoretically  and  experimentally  (see  e.g.  the  sur¬ 
vey  paper  by  Talbot,  Ref.  3l)>  it  can  be  assumed  that  the  one-dimensional  shock 
structure  in  a  pure  gas  is  reasonably  well  described  (at  least  experimentally). 
In  connection  with  early  experimental  shock  structure  studies,  the  Fh.D.  work 
of  F.S.  Sherman  (Ref.  32)  deserves  some  mentioning  at  this  point.  In  his  re¬ 
port  (and  thesis)  he  describes  the  development  and  operation  of  various  shock 
holders.  The  ideas,  which  he. put  forth  on  this  subject,  were  used  in  the  de¬ 
sign  of  the  shock  holder  described  in  Section  4.4.  Sherman  determined  tempera¬ 
ture  profilesthrough  steady  flow  shock  waves  (produced  by  shock  holders  placed 
in  a  supersonic  wind  tunnel  flow)  with  a  free-molecule  equilibrium  temperature 
probe.  Most  experimental  shock  structure  work  since  then  has  been  performed 
in  shock  tubes  and  shock  tunnels. 

For  the  present  work  the  normal  shock  wave  in  a  pure  gas  appear¬ 
ed  to  be  a  convenient  flow  configuration  for  testing  the  accuracy  of  the  elec¬ 
tron  beam  probe  under  the  extreme  conditions  of  strong  streamwise  gradients. 
There  existed  some  doubt,  whether  the  spatial  resolution  of  the  electron  beam 
probe  was  sufficient  to  give  a  true  density  profile  through  the  shock.  In 
this  context  it  should  be  kept  in  mind  that  the  spatial  resolution  in  the  flow 
direction  was  determined  by  the  electron  beam  width,  which  was  between  §mm 
and  1  mm.  Furthermore,  it  could  be  expected  that  the  distortion  of  the 
fluorescent  halo  (see  Section  5*3) >  caused  by  the  severe  density  gradient  in 
the  shock  resulted  in  a  possible  error  in  the  measured  density. 

To  obtain  an  idea  of  the  limits  imposed  by  the  spatial  resolu¬ 
tion  (or  the  lack  of  it)-  of  the  electron  beam  probe,  density  measurements  were 
made  through  plane  shock  waves  produced  by  means  of  the  shock  holder,  which 
was  placed  into  fre.e  jets  of  pure  argon  and  pure  helium.  Figure  5*l4  is  a  flow 
visualization  photograph  of  the  flow  configuration  obtained  by  mapping  out  the 
density  field  with  the  electron  beam.  It  depicts  a  Mach  4.5  shock  in  argon. 

The  density  profile  through  this  shock  (along  the  axis  of  symmetry),  as  mea¬ 
sured  with  the  optical  analyzer,  is  given  in  Fig.  5*15. 
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The  flow  in  the  free  jet  is  not  uniform;  it  has  streamwise 
gradients,  and  the  streamlines  diverge.  It  is,  therefore,  difficult  to  esta¬ 
blish  a  close  correspondence  of  the  measured  shock  profile  with  that  of  a  one¬ 
dimensional  shock  transition.  On  the  other  hand,  an  approximate  comparison 
can  be  made  by  referring  the  measured  density  profile  to  a  free  stream  base, 
which  follows  closely  the  axial  density  profile  of  the  undisturbed  free  jet 
(the  iseritropic  expansion  curve  for  the  free  jet  has  been  used  here) .  On  the 
downstream  side  of  the  shock  this  reference  curve  is  allowed  to  asymptotically 
approach  a  constant  downstream  value,  which  is  determined  in  such  a  way  as  to 
give  the  correct  theoretical  density  ratio  across  a  normal  shock  of  a  given 
Mach  number  (see  Fig.  5*15) •  An  approximate  shock  Mach  number  has  been  in¬ 
ferred  from  the  position  of  the  shock  in  the  free  jet  expansion,  which  has 
been  assumed  to  be  close  to  isentropic. 

Previous  experimental  shock  thickness  measurements  made  in 
argon  (see  e.g.  Ref.  33  and  3*0  indicate  that  the  one-dimensional  maximum  slope 
thickness,  based  on  the  density  profile,  ranges  between  4  and  5  upstream  mean 
free  path  lengths  for  shock  Mach  numbers  between  3  and  10.  For  the  argon  shock 
in  Figs.  5.14  and  5*15  the  equivalent  upstream  mean  free  path,  based  on  the 
asymptotic  value  of  the  reference  density  curve, is  about  1  mm.  On  this  basis, 
the  experimental  maximum  slope  thickness  (6S)  of  6  mm  ,  as  measured  with  the 
electron  beam,  appears  to  be  about  30$  too  large.  Although  this  discrepancy 
is  still  within  the  scatter  of  the  experimental  shock  thicknesses  reported  in 
the  literature,  it  is  perhaps  an  indication  that  the  spatial  resolution  of 
the  electron  beam  densitometer  prohibited  the  exact  analysis  of  shocks,  which 
were  thinner  than  about  1  cm.  On  the  other  hand,  this  discrepancy  could  just 
as  well  be  explained  by  the  fact  that,  strictly  speaking,  the  shock  was  not 
one-dimensional,  because  of  shock  curvature  with  respect  to  the  streamlines, 
non-uniform  upstream  flow  conditions,  and  the  finite  diameter  of  the  shock 
holder.  From  this  point  of  view  there  is  thus  no  conclusive  evidence  that 
the  limits,  imposed  by  the  spatial  resolution  of  the  electron  beam  probe,  re¬ 
sult  in  any  artificial  broadening  of  the  shock  transition.  Since  all  other 
shock  waves  discussed  in  this  report  are  wider  than  6  mm,  this  matter  has 
been  given  no  further  consideration. 

The  measured  density  profile  through  a  Mach  8.5  shock  in  helium 
is  given  in  Fig.  5«l6.  In  this  case  the  "equivalent"  maximum  slope  thickness 
of  the  shock  of  19  mm  (as  compared  to  an  upstream  mean  free  path  of  4.7  mm) 
is  in  good  agreement  with  other  results  reported  in  the  recent  literature 
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5.6  Static  Calibration  of  Gas  Mixtures  with  the  Aid  of  a  Mass-Spectrometer 

The  fact  that  secondary  processes  are  responsible  for  spread¬ 
ing  the  light  emission  in  a  pure  gas,  strongly  suggests  the  possibility  that 
the  helium  signal  observed  in  a  gas  mixture  might  be  influenced  by  the  pre¬ 
sence  of  argon  and  vice  versa.  This  could  come  about  by  the  transfer  of  pri¬ 
mary  excitation  energy  from  one  atomic  species  to  the  other  by  radiative  reson¬ 
ance  coupling,  by  interatomic  collisions,  or  by  secondary  electrons.  Radiative 
resonance  coupling  can  be  eliminated  immediately  on  the  grounds  that  the  energy 
levels  of  the  excited  states  of  helium  are  widely  separated  from  those  of 
argon.  It  turns  out,  that  all  excited  Arl  states  are  more  than  4  ev  below  the 
first  excited  state  of  helium.  Furthermore,  all  excited  Aril  states  are  at 
least  4.5  ev  above  the  ionization  limit  of  Hel  and  more  than  20  ev  below  the 
first  excited  state  of  Hell.  This  state  of  affairs  makes  an  inter- species 
energy  transfer  ,y  atomic  collisions  also  very  unlikely,  except,  perhaps,  by 
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some  rare  ’’collisions  of  the  third  kind",  in  which  the  colliding  atoms  are 
both  in  excited  states.  This  leaves  the  possibility  of  cross-excitation  by 
secondary  electrons . 

A  series  of  experiments,  using  gas  mixtures  of  known  composition 
and  at  known  densities,  were  performed  to  test  for  any  cross-excitation.  The 
helium-argon  mixtures  were  again  introduced  into  the  tunnal  at  point  A  (Fig*  4.1). 
Separate  helium  and  argon  supplies  were  used.  The  flow  rate  of  each  gas  could 
be  accurately  controlled  and  was  measured  by  a  separate  flow  rater.  Since  at 
low  tunnel  pressures  (below  50  p  Hg)  the  diffusion  pumps  pump  away  the  helium 
faster  than  the  argon,  it  follows  that  the  helium  to  argon  ratio  in  the  tunnel 
was  not  identical  with  the  ratio  of  the  flow  rates  into  the  tunnel.  It  was, 
therefore,  necessary  to  obtain  an  independent  measurement  of  the  gas  composi¬ 
tion  in  the  tunnel.  This  was  done  by  sampling  the  stagnant  gas  mixture  from 
a  point  near  the  electron  beam  with  a  free  molecule  extraction  tube  (0.15  mm 
dia.  x  5  cm  lg.  hypodermic  stainless  steel  tubing).  The  sampling  tube  was 
connected  to  a  "VEEC0-MS-9  AB"  mass  spectrometer  type  leak  detector.  The 
experimental  arrangement  is  shown  schematically  in  Fig.  5.17.  In  all  tests 
performed  with  this  system  the  pressure  behind  the  sampling  tube,  as  measured 
by  a  hot  cathode  ionization  gauge,  was  relatively  independent  of  the  tunnel 
pressure  and  ranged  between  5  x  10“ 5  torr  and  10_4  torr  at  all  times.  The 
helium  ion  current  obtained  from  the  mass  spectrometer  was  calibrated  against 
helium  pressure  in  the  test  section  by  admitting  pure  helium  into  the  tunnel 
and  by  measuring  the  helium  pressure  accurately  with  a  McLeod  gauge.  This 
calibration  curve  is  reproduced  in  Fig.  5*18. 

The  calibrated  system  was  then  used  to  determine  the  partial 
pressure  of  helium  in  the  gas  mixtures  studied.  The  calibration  curve  for  the 
helium  ion  current  could  be  used  directly  to  this  effect,  provided  the  flow 
through  the  sampling  tube  was  free  molecular,  in  which  case  the  flow  of  argor 
atoms  through  the  tube  did  not  influence  the  flow  of  helium  atoms  through  the 
tube.  For  the  case  of  helium  the  Knudsen  number  based  on  the  mean  free  path 
in  the  tunnel  and  the  sampling  tube  diameter  was  larger  than  5  for  tunnel 
pressures  below  200  (i  Hg;  for  the  case  of  argon  this  was  so  for  pressures 
below  110  p  Hg/  The  molecular  mean  free  path  in  the  1"  copper  tube*  connect¬ 
ing  the  sampling  tube  to  the  leak  detector  was  larger  than  60  cm  at  all  times. 

The  mixture  composition  in  the  tunnel  was  varied  from  1$  helium 
to  90$  helium,  and  the  total  pressure  in  the  test  section,  which  was  measured 
with  the  McLeod  gauge,  ranged  from  5  P  Hg  to  150  p  Hg.  The  helium  and  argon 
signals  from  the  electron  beam  densitometer  were  calibrated  for  these  helium- 
argon  mixtures  as  follows:  Through  a  control  valve  and  a  flow  rater  pure 
helium  was  admitted  into  the  tunnel  while  one  or  two  oil  booster  pumps  kept 
the  tunnel  pressure  at  an  equilibrium  value.  The  test  section  pressure 
(McLeod),  the  helium  ion  current  from  the  mass  spectrometer,  and  the  helium 
and  argon  signals  from  the  optical  analyzer  were  recorded.  Then  a  specified 
amount  of  argon  was  added  to  the  inflow  through  a  second  control  valve.  The 
helium  flow  rate  was  adjusted  until  the  helium  ion  current  from  the  mass 
spectrometer  was  again  at  the  previously  recorded  value  (i.e.  the  partial 


*  This  tube  had  been  thoroughly  cleaned  with  hot  detergent,  with  acetone, 
and  with  chlorothene,  and  had  been  baked  out  with  a  heat  gun  at  about  300°C 
to  reduce  outgassing. 
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pressure  of  helium  in  the  tunnel  was  the  same).  The  new  test  section  pressure 
and  the  electron  beam  densitometer  readings  were  recorded.  Then  the  inflow  of 
argon  was  increased  again  following  the  same  procedure.  In  this  manner  the 
partial  pressure  of  argon  was  increased  step  by  step,  while  keeping  the 
partial  pressure  of  helium  constant,  until  the  total  pressure  exceeded,,  lOOnHg. 
At  the  end  of  each  run  the  pure  helium  calibrations  were  checked  again  at  the 
particular  helium  pressure. 

The  procedure  described  above  offered  the  most  direct  way  of 
testing  for  any  influence  of  the  argon  on  the  helium  radiation  measured  by 
the  optical  analyzer.  Several  runs  were  made  starting  with  different  helium 
pressures.  After  subtracting  the  photomultiplier  dark  currents  and  the  error 
signals  it  was  found  that  within  the  experimental  error  the  helium  radiation 
at  5016a  remained  constant  for  each  run,  independent  of  the  amount  of  argon 
present . 


In  order  to  test  whether  the  presence  of  helium  had  any  in¬ 
fluence  on  the  argon  fluorescence,  the  same  data  could  be  used.  For  each  test 
the  partial  pressure  of  argon  in  the  test  section  was  obtained  by  subtracting 
the  partial  pressure  of  helium  (leak  detector)  from  the  total  pressure  (McLeod) . 
By  comparing  the  argon  signals  obtained  in  the  helium-argon  mixtures  with  the 
corresponding  pure  argon  calibration  values,  no  influence  of  the  helium  on  the 
argon  signal  could  be  detected. 

The  experimental  results  of  these  tests  are  summarized  in  Figs. 
5.19  and  5.20.  The  magnitude  of  the  probable  error  is  indicated  by  the  ver¬ 
tical  "error  bar"  drawn  through  each  experimental  point.  These  errors  are 
based  on  an  estimate  of  the  following  effects,  which  are  believed  to  be  re¬ 
sponsible  for  most  of  the  scatter  of  points  in  Figs.  5*19  and  5*2C: 

(1)  Small  fluctuations  of  the  tunnel  pressure  occurred  during  each  test 
because  of  the  somewhat  uneven  pumping  of  the  oil  vapour  booster  pumps  on 
the  tunnel.  Owing  to  their  different  time  constants  the  instruments  used  in¬ 
dicated  different  average  densities  or  pressures.  The  time  constants  involved- 
were: 

electron  beam  densitometer  t  <  1  sec 
McLeod  gauge  t  «  1  min 

mass  spectrometer  3  min<T  <  6  min 

(2)  The  ion  current  amplifier  in  the  leak  detector  showed  signs  of  drift 
and  of  instability  with  respect  to  power  line  fluctuations. 

It  is  interesting  to  note  that,  in  order  to  keep  the  helium 
pressure  constant  during  a  run,  the  helium  inflow  into  the  tunnel  had  to  be 
reduced  as  more  argon  was  added.  This  can  be  attributed  to  the  fact  that  above 
30  p  Hg  the  volumetric  pumping  speed  of  the  booster  pumps  decreased  with  in¬ 
creasing  pressure.  At  lower  pressures  this  must  be  attributed  to  a  decrease 
in  the  preferential  pumping  of  helium  with  increasing  pressure.  The  latter 
effect  has  been  investigated  by  leaking  a  constant  composition  gas  mixture 
(57*5$  He  -  42.5 1o  Ar)  from  a  high  pressure  cylinder  into  the  tunnel  and  by 
measuring  the  equilibrium  composition  in  the  tunnel  with  the  mass-spectrometer 
for  various  tunnel  pressures.  The  results  are  plotted  in  Fig.  5»21,  which 
seems  to  indicate  that  the  effect  has  virtually  disappeared  at  pressures  above 
50  p  Hg. 
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Since  the  tests  discussed  in  this  section  indicate  r.o  noticeable 
coupling  between  the  helium  radiation  and  the  argon  radiation,  it  must  be  con¬ 
cluded  that  the  secondary  electrons ,  which  originate  in  the  electron  beam 
from  ionizing  collisions  with  one  atomic  species,  have  on  the  average,  not 
enough  kinetic  energy  to  excite  atoms  of  the  other  species.  In  this  connection 
it  .is  of  interest  to  note  that  for  17 .5  kv  electrons  the  ionization  cross- 
section  of  the  argen  atom  is  about  ten  times  that  of  the  helium  atom  (see  Ref. 
35) »  In  mixtures  containing  more  than  10%  argon,  therefore,  the  majority  of  the 
secondary  electrons  produced  result  from  argon  ionization  events.  Hence,  the 
data  presented  in  Figs.  5.19  and  5-20  provide  strong  evidence  for  the  conclu¬ 
sion  that  secondary  electrons  originating  from  argon  atoms  have  energies  of 
less  than  20  ev.*  On  the  other  hand,  they  provide  relatively  weak  evidence 
for  the  secondary  electrons  originating  from  helium  atoms  to  have  er.ergr.es 
of  less  than  12  ev  .** 


VI. 


DISCUSSION  OF  SYSTEMATIC  ERRORS 


On  the  basis  of  the  preliminary  calibration  and  validity  tests 
the  following  remarks  can  be  made  about  some  of  the  possible  sources  of  error 
in  the  final  results. 


(l)  If  r.o  spurious  background,  like  is  present,  the  helium  and  argon 
"error  signals"  due  to  cross-talk  between  the  optical  filters  present  no  pro¬ 
blem,  since  these  can  be  accurately  calibrated  and  subtracted  from  the  final 
readings.  Owing  to  their  relative  smallness  (see  Fig.  5°2),  these  error  sig¬ 
nals  become  significant  only  for  the  rare  component  in  mixtures  containing 
very  small  amounts  of  either  hc-lium  or  argon. 


(2)  In  a  high  speed  flow  the  downstream  drift  of  metastatle  argon  atoms , 
resulting  in  a  partial  "less"  of  the  argon  halo,  causes  the  argon  signal  to 
be  lower  than  the  calibrated  value  for  a  given  density.  The  free  jet  tests 
in  pure  gases  show  this  effect  to  be  small,  and  it  has  therefore,  been  ignored. 


(3)  Simple  density  gradients  across  or  along  the  electron  beam  "distort" 
the  halo.  The  observed  signal  is,  however,  not  appreciably  afferted,  since 
the  effects  of  the  low  density  region  on  the  one  side  approximately  compen¬ 
sate  for  the  effects  of  the  high  density  region  on  the  other  side,  Second  spa¬ 
tial  density  derivatives  across  or  along  the  electron  beam  may,  on  the  other 
hand,  influence  the  observed  signals  significantly.  Second  density  derivatives 
across  the  beam  are  not  serious,  provided  almost  all  of  the  light  in  tne  halo 
is  collected  for  the  analysis.  Those  occurring  along  the  beam,  however,  in¬ 
variably  introduce  an  error  in  the  measured  density.  In  the  present  case  the 
most  severely  affected  measurements  are  those  made  on  the  axis  of  the  f’ec  jets 
close  to  the  5  mm  orifice.  Here  a  high  density  region  is  closely  suv rounded 
by  a  region  of  low  density,  resulting  in  a  diminished  halo  intensity.  For 
the  argon  density  measurements  small  corrections  have  been  applied  in  these 


*  Energy  above  the  ground  state  of  the  first  excited  state  of  the  helium 
atom  is  19.86  ev. 

**  Energy  above  the  ground  state  of  the  first  excited  state  of  the  argon 
atom  is  11,5  ev. 


27 


severe  cases,  based  on  an  estimate  of  the  effect.  For  the  helium  this  effect 
has  been  ignored. 

(4)  The  mass-spectrometer  tests  have  shown  that  any  coupling  between  the 
helium  and  argon  fluorescence,  through  secondary  electron  excitation  (or  any 
other  mechanism),  does  not  introduce  any  significant  error  into  the  partial 
density  measurements  of  the  component  gases  in  helium-argon  mixtures. 


(5)  The  reflection  and  interception  of  light  by  the  nozzle  face  and  the 
shock  holder  surfaces  has  been  found  to  affect  the  density  readings  near 
these  objects  (see  Section  5»l)*  This  effect  has  been  calibrated  for  the 
relevant  spectral  lines  in  the  density  range  encountered.  Corresponding 
correction  factors  have  been  applied. 

(6)  The  dark  currents  of  the  photomultiplier  tubes  were  measured  before 
and  after  every  run  and  have  been  subtracted  from  every  reading.  The  magni¬ 
tudes  of  the  average  dark  currents  were 

5  x  10"^®  amp  for  the  EMI- 9502  S  tube,  and 

2  x  10" 9  amp  for  the  EMI-9558A  tube  . 

In  most  instances  these  currents  were  completely  negligible  when  compared  to 
the  normal  signal  strengths  (see  Fig.  5.l)* 

At  this  point  a  word  is  in  order  about  the  possibility  of  con¬ 
densation  occurring  in  the  free  jets,  since  this  might  complicate  the  free 
jet  structure  and  influence  the  diffusion  processes  taking  place  in  the  free 
jet  and  the  shock  waves  produced  therein.  In  the.  free  jet  expansions  (stag¬ 
nation  temperature  was  about  300°K)  the  local  temperatures  and  pressures  de¬ 
creased  such  that  under  equilibrium  conditions  the  argon  should  have  conden¬ 
sed  (into  the  solid  state)  between  Mach  numbers  of  4  and  5  in  the  jet.  Owing 
to  the  low  densities  and  the  high  rates  of  expansion  occurring  in  these  free 
jets,  however,  condensation  does  not  easily  take  place,  and  a  high  degree  of 
supersaturation  can  be  obtained.  For  nucleation  to  occur  in  the  free  jet 
"Arg-molecules"  have  to  be  formed  first  by  triple  collisions  (a  third  particle 
is  needed  to  carry  away  the  heat  of  sublimation).  Anderson,  Andres,  Fenn, 
and  Maise  (Ref.  36)  have  made  a  calculation,  based  on  kinetic  theory,  of  the 
mole  fraction  of  Ar2  to  be  found  in  free  argon  jets  expanding  from  room  tem¬ 
perature  to  the  collisionless  state.  The  predicted  values  are  in  agreement 
with  experimental  values  obtained  by  Milne  (Refs.  37  and  38),  who  observed 
"cluster  populations"  in  free  jets  using  a  mass  spectrometer.  On  the  basis 
of  these  studies  it  is  estimated  that  the  mole  fraction  of  Ar2  in  any  of  the 
jets,  used  in  the  present  work,  is  below  5  x  10”^.  Accordingly,  it  is  believed 
that  condensation  effects  can  be  ignored. 
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VII .  FREE  JET  RESULTS  (Helium-Argon  Mixtures) 


Having  gained  sufficient  confidence  in  the  technique  by  the 
foregoing  tests  the  optical  electron  beam  probe  was  used  to  study  the  density 
distribution  in  the  free  jet  expansion  of  helium-argon  mixtures.  The  mole 
fraction  of  argon  in  the  stagnation  chamber  ranged  from  0.05  to  0.5*  The  gas 
mixtures  were  bought  in  pressurized  cylinders  from  the  Union  Carbide  Corpora¬ 
tion.  Although  the  gas  composition  in  these  cylinders  was  known  approximately, 
a  separate  quantitative  analysis  (see  Appendix  A)  of  these  gas  mixtures  had 
to  be  performed,  since  the  reduction  of  the  free  jet  data  required  an  accurate 
knowledge  of  the  f0-values.  The  cylinders  used  were  found  to  contain  the 
following  mole  fractions  of  argon. 


cylinder  1 
cylinder  2 
cylinder  3 
cylinder  4 
cylinder  5 


*o 

**0 

*o 

*0 


0.12 

0.21 

0.425 

0.11 

0.49 


0.005 


Other  mixtures  were  metered  into  the  tunnel  from  separate  gas  supplies.  The 
stagnation  temperature  was  approximately  300°K  for  all  runs.  The  jet  Reynolds 


numbers  varied  between  100  and  10,000. 
as  follows: 

D  =  47.5  mm 
1)  =  15  mm 
D  =  5  mm 


The  range  covered  by  each  nozzle  was 

100  <  Re0  <  500, 

200  <  Re0  <  2000, 

500  <  Re0  <  10,000. 


The  Reynolds  number  ranges  overlap  sufficiently,  so  that,  for  some  free  jet 
runs,  Reynolds  numbers  could  be  duplicated  with  different  orifices. 


It  may  be  argued,  on  account  of  the  non-zero  value  of  0  -fQ 
at  the  center  of  the  orifice  (Fig.  2.2),  that  the  mole  fractions  of  the  gas 
components  in  the  stagnation  .chamber  could  well  be  different  from  those  in 
the  supply  cylinder.  This  would  indeed  be  true  if  the  integral  of  0-fo  , 
integrated  over  the  orifice  area,  were  non- zero.  For  the  reduction  of  the 
experimental  data  it  has  been  assumed,  that  in  the  orifice  the  value  of 
0-f 0  is  sufficiently  negative  off  the  axis  to  exactly  balance  the  positive 
contribution  on  the  jet  axis.  An  argument  is  presented  in  Appendix  B  to 
support  this  assumption.  Consequently,  the  stagnation  chamber  value  of  f0 
has  been  taken  to  be  identical  with  that  in  the  supply  cylinder  (or  with  the 
value  derived  from  the  volume  flow  rates  metered  in) . 


Figure  5»13  shows  the  general  structure  of  the  jet  flow  fields. 
As  is  evident  from  these  photographs,  the  lower  Reynolds  number  jets  exhibit 
a  very  diffuse  shock  envelope.  For  these  jets  the  barrel  shock  and  the  Mach 
disc  have  become  broadened  to  the  point,  where  the  extent  of  the  undisturbed 
central  flow  field  is  considerably  reduced.  In  agreement  with  the  photo¬ 
graphs  the  axial  density  surveys, performed  with  the  electron  beam  probe, 
indicate  a  density  deviation  from  the  isentropic  behaviour  somewhat  ahead  of 
the  expected  Mach  disc  location. 


An  attempt  has  been  made  to  account  for  the  different  appear¬ 
ances  of  these  jets  in  terms  of  the  flow  parameters.  The  general  features 
of  an  underexpanded  jet  are  characterized  by  the  size  of  the  first  expansion 
cell  and  by  the  diffuseness  of  the  barrel  shock  and  Mach  disc.  The  distance 
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,  .w.  ±u<=  itAdfcive  tmcK ness  of  the  barrel  shock  at  any  point, 

or  that  of  the  Mach  disc,  may  be  defined  by  the  ratio  -formed  by  the  shock 
thickness,  5  ,  and  the  streamline  distance,  s,  from  the  orifice  plane  to  the 
center  of  the  shock.  If  it  is  assumed  that  the  density  along  the  jet  axis 
falls  off  as  (x/3))'8  (corresponding  to  a  simple  source  flow  model  with  the 
source  located  near  the  orifice),  and.  if  the  shock  thickness  is  assumed  to  be 
equal  to  4.5  times  the  distance  of  a  molecular  mean  free  path  ahead  of  the 
shock  (a  value,  which  agrees  reasonably  well  with  experimental  data  for  M+>  4), 
then  the  relative  thickness  of  the  Mach  disc  can  be  shown  (see  Appendix  C]  to 
vary  inversely  as  the  jet  Beynolds  number  and  directly  as  the  jet  pressure 
ratio  raised  to  a  power,  which  falls  between  l/3  and  l/2,  depending  on  the 
viscosity-temperature  law  assumed.  If  -  (Tj/T0)2°  is  assumed,  then  the 

relative  thickness  of  the  Mach  disc  becomes, 


(7+1)2  (7-1) 


(13) 


This  equation  applies  to  a  pure  gas.  For  a  gas  mixture  the 
shock  may  be  several  times  wider.  When  the  jet  Beynolds  number  and  the  jet 
pressure  ratio  are  such  that  5/xj,j  becomes  approximately  unity,  then  no 
region  in  the  jet  remains  undisturbed  by  the  ambientxgas.  For  the  barrel 
shock  6/s  can  be  considered  proportional  to  (pc/Poo)^  Re^  in  a  very  approxi¬ 
mate  sense  onlyc  A  qualitative  analysis  of  a  number  of  flow  visualization 
photographs  confirms  the  above  dependence. 

Some  of  the  experimental  density  profiles  obtained  along  the  axes 
of  various  free  jets  are  presented  in  Figs.  7*1  to  7.8.  In  these  figures  the 
partial  number  densities  of  helium  and  argon  atoms,  normalized  with  respect 
to  the  corresponding  partial  number  densities  in  the  stagnation  chamber,  have 
been  plotted  against  axial  distance  in  terms  of  x/D.  The  original  mole  frac¬ 
tion  of  argon,  the  jet  Reynolds  number,  the  orifice  diameter,  and  the  pressures 
on  either  side  of  the  orifice,  are  stated  for  each  run.  The  gas  mixtures  as  a 
whole  follow  the  Isentropic  expansion  curve  quite  well,  which  drops  off 
approximately  as  the  inverse  square  of  the  axial  distance  for  x  larger  than 
2D.  The  first  experimental  downstream  deviation  from  this  dependence  for 
any  one  run  indicates  the  point  at  which  the  surroundings  start  to  affect  the 
axial  flow.  The  apparent  discrepancy  between  the  shape  of  the  isen trope  and 
the  experimental  curves,  very  close  to  the  orifice  (x/D  <  0.5),  is  simply  a 
consequence  of  the  fact,  that  the  isentropic  expansion  through  a  sonic  orifice 
is  not  exactly  known.  In  all  cases  the  argon  curve  is  above  the  helium  curve, 
indicating  an  argon  enrichment  along  the  jet  centerline.  It  is  evident  from 
these  figures  that  the  axial  argon  enrichment  decreases  with  increasing  jet 
Reynolds  number.  As  shown  by  Figs.  7*7  end  7.8  the  diffusive  component  separa¬ 
tion  on  the  jet  axis  becomes  very  small  indeed  for  Reynolds  numbers  of  2000 
and  3000.  For  jet  Reynolds  numbers  above  5000  the  separation  was  smaller  than 
the  experimental  error  and  could  not  be  detected.  For  all  density  profiles 
shown,  the  density  increase  caused  by  the  "Mach  disc”  (i.e.  the  diffuse  shock 
envelope  closing  in)  is  clearly  evident.  As  predicted  by  one-dimensional  shock 
theory  the  helium  compresses  ahead  of  the  argon,  resulting  in  a  helium  enrich¬ 
ment  in  the  Mach  disc.  The  corrected  and  smoothed  out  experimental  density 
data  are  believed  to  be  accurate  to  within  one  or  two  percent. 
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From  these  experimental  density  profiles,  the  local  mole  frac¬ 
tions  of  argon  can  easily  be  determined.  For  the  free  jets  starting  out  with 
original  mole  fractions  ofr argon  of  0.12  and  0.425  these  mole  fractions  nave 
been  plotted  against  (x/D)e  for  various  Reynolds  numbers  in  Figs.  7»9  and 
7.10.  The  solid  curves  give  the  theoretical  mole  fractions  as  per  Fig.  2.2 
The  downstream  "tails",  indicating  a  sudden  decrease  in  the  mole  fraction  of 
argon,  are  caused  by  the  Mach  discs.  It  should  be  noted  that  free  jets 
originating  from  different  stagnation  conditions  and  emerging  through  differ¬ 
ent  sized  orifices,  but  having  approximately  the  same  Reynolds  number  (and  the 
same  fQ),  show  identical  species  separation  effects.  Also,  in  the  undisturb¬ 
ed  portions  of  the  free  jet  expansions  the  agreement  with  theory  is  quite 
reasonable. 


In  Fig.  7*11  the  experimental  data  for  gas  mixtures  ranging 
from  5$  to  argon  and  for  Reynolds  numbers  between  100  and  1000  have  been 
summarized  and  are  compared  with  Sherman's  theory  (Ref.  8).  Only  those  por¬ 
tions  of  the  free  jets  which  were  not  affected  by  the  shock  envelope  have 
been  used  in  this  graph.  Note  that  in  this  representation  the  experimental 
errors  are  multiplied  by  the  Reynolds  numbers  of  the  jets,  resulting  in  a 
fair  amount  of  experimental  scatter.  It  is  evident  that  Sherman's  first 
order  approximation  predicts  the  actual  species  separation  on  the  jet  axis 
cuite  well,  even  for  Reynolds  numbers  as  small  as  100. 

Point-by-point  measurements  of  the  partial  densities  of  helium 
and  argon  were  made  throughout  the  free  jet,  off  the  axis  as  well  as  on  the 
axis,  for  two  cases: 

Jet  1:  fQ  =  0.12,  Re0  =  533,  D  =  15  mm,  P0/Poo  =  150  . 

Jet  2:  fQ  =  0.21,  Re0  =  224,  D  =  15  mm,  Pq/p^  =  188  . 

For  jet  1  the  radial  number  density  profiles  at  various  axial  positions  are 

given  in  Fig.  7*12.  The  mole  fraction  of  argon  throughout  the  jet  is  shown 

in  Fig.  7.13*  The  "depression"  in  the  downstream  direction  after  the  initial 
increase  in  f  is  caused  by  the  Mach  disc.  A  pronounced  argon  deficiency 
exists  also  near  the  edge  of  the  jet.  This  is  believed  to  be  partially  due 
to  the  presence  of  the  barrel  shock  and  partially  due  to  the  fact  that  radial 
pressure  diffusion  causes  an  argon  enrichment  on  the  axis  and  a  corresponding 
argon  deficiency  off  the  axis  of  the  jet.  The  gas  in  the  test  section  outside 
of  the  jet  is  slightly  more enriched  in  argon  than  the  stagnation  chamber  gas, 
because  the  diffusion  pumps  remove  the  helium  at  a  somewhat  faster  rate.  A 
flow  visualization  of  jet  1  is  given  by  the  bottom  left  photograph  of  Fig. 

5.13.  The  radial  partial  density  profiles  for  jet  2  are  presented  in  Fig.  7«l4; 
the  mole  fractions  of  argon  throughout  the  flow  field  are  plotted  in  Fig.  7.15. 
The  results  indicate  the  same  general  features  for  jet  2  as  were  observed  for 
jet  1.  Throughout  the  jet  the  separation  effects  are  somewhat  more  pronounced 
because  of  the  lower  Reynolds  number.  Furthermore,  the  features  which  can  be 
related  to  the  barrel  shock  and  Mach  disc  are  considerably  more  diffuse. 
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vm.  SHOCK  WAVE  RESULTS  (Helium-Argon  Mixtures) 

Normal  shock  waves  were  produced  in  front  of  the  partially 
choked  shock  sholder  p3.aced  in  the  free  jets.  Figure  8.1  shows  the  general 
appearance  of  the  shock  waves  investigated.  It  is  a  flow  visualization  photo¬ 
graph  of  a. Mach  9  shock  in  a  49 %  argon  mixture.  The  shock  appears  to  be 
quite  plane,  perhaps  slightly  "dished  in".  It  should  be  kept  in  mind,  however, 
that  a  plane  shock  in  a  diverging  flow  corresponds  to  a  curved  convex  shock  in 
a  parallel  flow.  From  this  point  of  view  the  shock  in  Fig.  8.1  is  not  dished 
in  enough  to  be  perpendicular  to  the  incoming  streamlines  at  all  points  and 
corresponds,  therefore,  to  a  shock  which  has  a  small  convex  curvature  with 
respect  to  the  streamlines.  This  is  also  true  for  the  other  shock  waves 
studied.  Furthermore,  the  shock  diameter  is  only  about  three  times  the  shock 
thickness.  The  shocks  can,  thus,  hardly  be  considered  one-dimensional. 

Figures  8.2  to  8.4  give  the  normalized  partial  density  profiles 
along  the  axes  of  the  underexpanded  jets  and  through  the  shock  waves  for  49 %, 

11%  and  1%  original  argon  mixtures,  respectively.  The  corresponding  profiles 
of  "f"  are  given  in  Figs.  8.5  to  8.7.  The  initial  argon  enrichment  is  caused 
by  the  free  jet  expansion.  In  all  eases  the  helium  compresses  ahead  of  the 
argon  in  concord  with  theory.  Radial  diffusion  in  and  behind  the  shock  is 
believed  to  be  responsible  for  the  pronounced  argon  surplus  inside  the  shock 
holder.  For  the  1%  argon  case  the  general  features  are  as  for  the  other 
cases,  except  for  the  behaviour  of  the  argon  density  in  the  upstream  portion 
of  the  shock.  The  partial  argon  density,  which  originally  is  significantly 
above  the  isentropic  expansion  curve,  decreases  below  the  isentrope  prior  to 
compression.  This  can  be  interpreted  as  an  acceleration  of  the  argon  atoms 
beyond  the  free  stream  value,  in  which  case  it  would  verify  the  effect  des¬ 
cribed  by  Sherman  (Ref.  9).  It  should  be  noted,  however,  that  the  scatter 
of  the  experimental  argon  densities  in  Fig.  8.4  is  relatively  large  because  of 
the  smallness  of  the  argon  signals  measured  in  this  case.  Furthermore,  the 
argon  error  signals  subtracted,  due  to  the  presence  of  helium,  were  larger 
than  the  true  argon  signals.  On  account  of  this  there  exists  some  uncer¬ 
tainty  about  the  argon  density  curve  drawn  in  Fig.  8.4.  The  effect,  however, 
appears  to  be  real,  and  it  remains  a  matter  of  interpretation,  whether  it  is 
ascribed  to  a  one -dimensional  acceleration  of  the  argon  or  to  a  three- 
dimensional  effect,  hitherto  unexplained. 

In  order  to  get  a  rough  comparison  of  the  experimental  shock 
profiles  with  the  existing  one-dimensional  theories  an  attempt  has  been  made 
to  compensate  for  the  fact  that  the  shock  is  embedded  in  a  flow  with  consid¬ 
erable  streamwise  gradients.  To  this  effect  an  arbitrary  reference  curve 
has  been  drawn  in  each  of  Figs.  8.2,  8,3  and  8.4.  This  reference  curve 
follows  the  isentropic  expansion  curve  into  the  shock  and  then  levels  out 
smoothly  toward  a  constant  downstream  reference  density,  which  was  chosen  to 
be  equal  to  the  isentropic  free  jet  density  at  the  centre  of  the  helium  shock. 

Of  course,  this  procedure  has  to  be  viewed  with  serious  reservations  and  is 
offered  here  only  in  an  attempt  to  compare  the  salient  features  in  a  qualitative 
way. 


Figure  8.8  gives  the  total  mass  density  divided  by  the  local 
reference  density  through  the  shock.  The  total  mass  density  has  been  calcu¬ 
lated  from  the  density  curves  for  the  individual  species.  The  distance  is 
given  in  terms  of  the  upstream  mean  free  path  defined  by 
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where  the  subscript  1  refers  to  conditions  just  ahead  of  the  shock.  Equation 
(l4)  is  the  accepted  formula  for  the  mean  free  path  of  smooth  elastic  spheres. 
Even  though  this  "mean  free  path"  has  no  direct  microscopic  meaning  for  a  gas 
mixture,  it  serves  as  a  convenient  reference  length.  In  calculating  this 
reference  length  the  shock  Mach  number,  M^_,  was  estimated  from  the  position 
of  the  shock  in  the  free  jet.  The  free  stream  quantities,  p]_,  T^,  and  ui 
were  then  obtained  from  p0,  T0  and  M-j_,  using  the  isentropic  relations.  The 
pure  helium  and  pure  argon  viscosities  at  Ti  were  estimated  according  to  the 
Chapman-Enskog  theory,  assuming  an  atomic  interaction  potential  of  the  Lennard- 
Jones  type.  Wilke's  formula  (Ref.  l6)  was  then  used  to  calculate  for  ‘tlie 
gas  mixture,  using  the  pure  gas  viscosities.  For  the  purpose  of  displaying 
several  shock  profiles  in  one  figure,  the  zero  point  of  the  axial  length  scale 
has  been  relocated  at  the  position,  at  which  the  total  mass  density  has  reach¬ 
ed  half  the  theoretical  density  increase  for  the  given  shock  Mach  number.  To 
distinguish  this  shifted  length  scale  from  the  previously  defined  x-scale  (with 
the  zero  point  at  the  orifice),  a  new  symbol,  x* ,  has  been  introduced  in  Fig. 
8.8. 


The  increase  in  the  non-dimensional  shock  thickness  with  in¬ 
creasing  mole  fraction  of  argon,  from  pure  helium  (curve  l)  to  50%  argon  (curve 
4),  is  evident  from  the  shock  profiles  in  Fig.  8.8,  which  are  all  for  an  up¬ 
stream  Mach  number  of  about  9*  The  compression  ratio  across  the  experimental 
pure  helium  shock  agrees  quite  well  with  one-dimensional  shock  theory  and 
provides  some  "justification"  for  the  use  of  the  somewhat  arbitary  reference 
curve.  In  the  case  of  the  gas  mixtures  the  total  mass  density  on  the  down¬ 
stream  side  of  the  shock  turns  out  to  be  too  high,  because  there  is  too  much 
argon  due  to  radial  diffusion  in  the  "curved"  shock. 

It  can  be  expected  that  both  the  shock  Mach  number  and  the 
Reynolds  number,  based  on  the  free  stream  conditions  and  the  radius  of  shock 
curvature  (i.e.  Rer  =  p^u^  rg/p-^),  become  scaling  parameters  for  the  radial 
diffusion  effect  behind  curved  shocks.  A  third  parameter  is  needed  to  describe 
the  dependence  of  the  effect  on  the  free  stream  composition.  Since  no  appro¬ 
priate  function  has  as  yet  been  derived  for  the  curved  shock  problem,  it  is 
assumed  here  that  the  function  Ej_,  based  on  the  free  stream  composition, 
adequately  describes  the  composition  dependence  of  the  radial  diffusion  effect. 
In  parallel  with  Sherman's  free  jet  analysis  one  might  then  expect,  that, 
for  a  given  shock  Mach  number,  the  quantity  (fg-f^Rer/E-^  is  approximately 
constant  for  all  curved  shocks  (fg  is  the  mole  fraction  of  argon  behind  the 
normal  part  of  the  shock).  Furthermore,  the  Mach  number  dependence  of  the 
argon  enrichment  is  probably  insignificant  for  M^_  >  5* 

For  the  experimental  shock  waves  discussed  here,  it  is  diffi¬ 
cult  to  estimate,  what  the  equivalent  shock  curvature  would  be,  if  the  shocks 
were  embedded  in  a  parallel  flow.  An  approximate  transformation  from  the 
diverging  flow  of  the  free  jet  to  a  parallel  flow  can  be  obtained,  if  one 
imagines  the  incident  streamlines  rigidly  attached  to  the  shock  surface,  and 
if  one  then  rotates  these  streamlines  into  a  parallel  configuration,  thereby 
forcing  the  shock  surface  into  a  more  convex  shape  (the  angle  of  incidence  of 
each  streamline  remaining  constant  in  the  process).  On  this  basis  it  has  been 
estimated  from  flow  visualization  photographs,  that  the  equivalent  radii  of 
shock  curvature  ranged  between  10  cm  and  20  cm,  and  that  the  quantity 
(f2-fi)Rer/Ei  is  of  the  order  of  100  for  these  shocks. 
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It  is  also  difficult  to  define  a  meaningful  shock  thickness 
for  the  shocks  of  Fig.  8.8,  since  the  downstream  halves  of  the  shock  waves  2, 
3,  and  k  differ  radically  from  any  corresponding  one- dimensional  shocks.  The 
upstream  halves  of  the  shock  waves,  on  the  other  hand,  should  not  be  signifi¬ 
cantly  influenced  by  radial  diffusion  and  should  correspond  closely  to  those 
of  one-dimensional  shocks.  This  appears  to  be  a  reasonable  assumption,  be¬ 
cause  in  the  supersonic  part  of  the  shock  the  downstream  influences  are  expect 
ed  to  produce  no  significant  curvature  of  the  streamlines.  An  approximate 
"equivalent  ID- shock  thickness",  6S,  has  therefore  been  defined  by  taking  the 
slope  of  the  total  mass  density  profile  at  x'  =0  and  extending  it  down  to 
p/pREF  =  1  and  UP  ^o  the  theoretical  density  ratio.  Thus  6S  is  defined  by 


6 


s 


(p2/p1)th  1 

[d(p/pREF)/dx' ]xt_o 


(15) 


The  values  of  6S  and  §s/\  obtained  in  this  fashion  are  tabu¬ 
lated  in  Fig.  8.8.  These  values  indicate,  that  for  a  shock  of  =  9  and 

^  =  0.5  the  shock  thickness,  in  terms  of  A]_,is  about  three  times  that  in  a 
pure  gas  (f^=  0  or  l). 


IX.  CONCLUSIONS 


Based  on  the  reported  experiments  and  discussionsthe  following 
concluding  remarks  can  be  made. 

(1)  The  method  of  using  a  thin  electron  beam  and  a  set  of  interference 
filters  for  a  spectroscopic  study  of  the  composition  of  rarefied  gases  has 
been  shown  to  give  accurate  results  for  argon-helium  mixtures.  Other  mix¬ 
tures  may  presumably  be  analyzed  in  a  similar  manner,  provided  that  the 
excitation  processes  involved  are  reasonably  well  understood.  It  appears  to 
be  necessary  to  make  preliminary  tests,  similar  to  those  outlined  in  Section 
V,  for  each  new  gas  mixture  and  for  each  new  set  of  spectral  lines  selected 

in  order  to  interpret  the  beam  intensity  measurements  correctly.  The  electron 
beam  provides,  at  the  same  time,  a  convenient  means  of  mapping  out  photo¬ 
graphically  the  density  field  under  study. 

(2)  No  radiative  coupling  between  the  helium  and  the  argon  fluorescence 
has  been  detected  in  helium-argon  mixtures.  Any  interspecies  transfer  of 
excitation  energy  by  secondary  electrons  must  evidently  be  insignificant. 

(3)  Sherman's  theory  for  the  species  separation  along  the  centerline  of 
an  underexpanded  free  jet  has  been  substantiated  experimentally  for  a  binary 
mixture  of  monatomic  gases.  According  to  this  theory  and  the  experimental 
tests  the  species  separation  in  the  undisturbed  free  jet  is  smaller  than  has 
previously  been  observed  by  investigators,  such  as  Stern  and  Waterman  or  Chow, 
who  analyzed  the  gas  sampled  by  a  skimmer  or  sampling  tube  placed  in  the  jet. 
The  species  separation  is  inversely  proportional  to  the  jet  Reynolds  number 
and  takes  place  almost  entirely  over  the  first  three  diameters  downstream  of 
the  orifice.  For  Reynolds  numbers  above  5000  the  species  separation  becomes 
insignificantly  small. 
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(4)  The  salient  features  predicted  for  the  one-dimensional  shock  struc¬ 
ture  in  a  binary  non-reacting  gas  mixture  have  been  observed.  For  all  shock 
waves  encountered  the  light  gas  has  been  found  to  compress  significantly  ahead 

of  the  heavy  gas,  resulting  in  an  enrichment  of  the  light  atomic  species  through¬ 
out  most  of  the  shock.  In  terms  of  the  upstream  gas  properties  the  shock 
transition  for  a  gas  mixture  has  been  observed  to  be  up  to  three  times  as  wide 
as  for  a  pure  gas .  For  a  Mach  9  shock  in  a  1$  argon  mixture  the  argon  has 
been  observed  to  expand  beyond  the  free  stream  value  in  the  upstream  portion 
of  the  shock  transition.  If  this  is  interpreted  as  a  one-dimensional  effect 
it  means  an  argon  acceleration  as  predicted  by  Sherman.  It  is  felt,  however, 
that  for  a  definite  proof  of  the  existence  of  this  effect,  is  is  necessary  to 
repeat  the  experiment  in  a  truly  uniform  one-dimensional  high  Mach  number  flow, 
which  was  not  attainable  with  the  present  UTIAS  low  density  facility. 

(5)  The  shock  wave  results  indicate  that  shock  curvature  may  produce 
very  significant  radial  diffusion  effects,  resulting  in  an  enrichment  of  the 
heavy  gas  on  the  centerline  behind  the  shock.  This  effect  might  be  used  to 
explain,  at  least  partially,  the  unrealistically  large  species  separation 
observed  when  sampling  tubes  are  introduced  into  a  supersonic  flow  field. 

There  remain  a  number  of  interesting  gas  dynamical  problems 
for  the  study  of  which  the  optical  electron  beam  densitometer  is  ideally 
suited.  A  short  list  of  the  most  obvious  extensions  of  the  present  experi¬ 
ments  is  given  hereunder. 

(.1)  A  study  could  be  made  of  the  diffusion  effects  in  the  free  jet  expan¬ 
sion  and  the  shock  compression  of  a  mixture  of  diatomic  gases  (say,  N2  and  Hg) . 
This  is  of  particular  interest,  since  the  rotational  temperatures  of  the  com¬ 
ponent  species  could  also  be  obtained  using  the  electron  beam  technique  (Ref. 

39).  The  component  temperatures  are  not  expected  to  be  equal  at  a  point  in¬ 
side  a  shock  wave,  for  example. 

(2)  Along  similar  lines,  a  study  of  the  diffusive  separation  in  gas  mix¬ 
tures,  the  components  of  which  having  different  specific  heat  ratios  (7), 
should  be  of  interest.  This  problem  also  remains  to  be  solved  theoretically. 

(3)  A  systematic  study  should  be  made  of  the  effect  of  shock  curvature 
on  the  degree  of  heavy  species  enrichment  behind  the  curved  shock. 

(4)  In  the  present  study  an  attempt  was  made  to  verify  one-dimensional 
shock  structure  theory  for  binary  gas  mixtures.  The  present  experiments 
have  shown  that  this  cannot  be  done  accurately  in  a  free  jet  flow  field. 

These  studies  should  be  repeated  in  a  truly  uniform  flow. 

(5)  Species  concentration  measurements  should  be  of  value  in  countless 
applications  involving  reacting  gas  mixtures. 

(6)  A  detailed  spectroscopic  study  of  the  radiation  resulting  from 
electron  beam  excitation  in  various  gases  is  of  interest  in  itself,  since  it 
may  lead  to  a  better  understanding  of  atomic  physics  and  to  improvements  in 
the  electron  beam  technique. 
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(7)  It  is  self-evident  that  the  optical  electron  beam  density  probe  can 
be  applied  to  numerous  gas  dynamical  problems  in  pure  gases  (e.g.  for  a  study 
of  boundary  layers,  wakes,  turbulence,  shocks  etc.). 

For  electron  beam  studies  in  helium-argon  mixtures  or  in  pure 
argon  the  spatial  resolution  can  be  improved  by  choosing  spectral  lines,  which 
do  not  show  a  fluorescent  halo.  In  retrospect  of  the  present  study  it  appears 
to  be  advantageous  to  sacrifice  intensity  in  favour  of  spatial  resolution, 
thereby  avoiding  needless  corrections.  From  that  point  of  view  the  Aril  lines 
between  A  4300a  and  A4500A  appear  to  be  suitable.  In  the  case  of  helium 
the  fluorescent  halo  observed  at  A  5016a  appears  to  have  little  effect  on 
the  spatial  resolution  except  at  very  low  densities. 
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FIG.  2.  2  THEORETICAL  MACH  NUMBER  VARIATION  AND  DIFFUSIVE 
SEPARATION  EFFECTS  ALONG  THE  AXIS  OF  AN  UNDER¬ 
EXPANDED  FREE  JET 


FIG.  2.  3  THEORETICAL  DEPENDENCE  OF  THE  REYNOLDS  NUMBER 

AND  OF  THE  PARAMETER  Re0/EC  ON  THE  GAS  COMPOSITION 
FOR  ARGON- HELIUM  MIXTURES 
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FIG.  4.2  EXTERNAL  VIEW  OF  GUN  ASSEMBLY  AND  INSTRUMENTATION 
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FIG.  4.4  DEFLECTING  PLATE  ASSEMBLY  AND  INNER  GUN  CHAMBER 
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FIG.  4.3  SCHEMATIC  DIAGRAM  OF  PUMPING  CONNECTIONS  AND  PROTECTIVE 
DEVICES 


FIG.  4.6  SPECTRAL  SENSITIVITY  CURVES  FOR  THE  PHOTOMULTIPLIER- 
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FIG.  4.  8  MODIFIED  VERSION  OF  THE  OPTICAL  ANALYZER  AS  USED 
FOR  MOST  OF  THE  EXPERIMENTAL  TESTS 


SUPPORT  TO 
ELECTRON  GUN  TUBE 


ENLARGED  VIEW  OF 
ENTRANCE  SECTIONS  OF 
OPTICAL  LIGHT  GUIDES 

FIBER  DIAMETER  =  0.076  mm 


VOLTAGE  DIVIDER 
BOX 
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MODIFIED  VERSION  OF  THE  OPTICAL  ANALYZER  AS  USED 
FOR  MOST  OF  THE  EXPERIMENTAL  TESTS 
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FIG.  4. 9 


WIRING  DIAGRAM  FOR  PHOTOMULTIPLIER  TUBES 


FIG.  4.10  SHOCK  HOLDER 
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FIG.  5. 1  PURE  GAS  CALIBRATION  CURVES  OF  PHOTOMULTIPLIER  CURRENT 
AGAINST  GAS  DENSITY 


FIG.  5.4 


ENERGY  LEVEL  DIAGRAM  FOR  NEUTRAL  HELIUM 
(taken  from  Ref.  24) 
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FIG.  5.  5  ELECTRON  BEAM  EXCITED  SPECTRUM  OF  HELIUM. 

THE  1/2  mm  WIDE  ELECTRON  BEAM  IS  AT  RIGHT  ANGLES 
TO  THE  SPECTROGRAPH  SLIT 
ELECTRON  BEAM:  17.  5  kv  -  200// amp 

HELIUM  PRESSURE:  100  //  Hg 
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FIG.  5.7  TABLE  OF  Ar  1  -  LINES  FOR  TRANSITIONS  FROM  THE 
3p5  4p  AND  3p5  5p  LEVELS  TO  THE  3p5  4s  LEVELS 
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FIG.  5.8  ELECTRON  BEAM  EXCITED  SPECTRUM  OF ^ARGON  IN  THE 
SPECTRAL  RANGE  FROM  3900  A  TO  4900  A. 

THE  SPECTROGRAM  WAS  OBTAINED  IN  AN  UNDEREXPANDED 
ARGON  JET  WITH  THE  I  mm  WIDE  ELECTRON  BEAM  AT 
RIGHT  ANGLES  TO  THE  SPECTROGRAPH  SLIT. 

THE  JET  FLOW  IS  FROM  LEFT  TO  RIGHT 

-  Spectrogram  by  D.  Marsden,  UTIAS  - 
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FIG.  5.  9  B  ELECTRON  BEAM  EXCITED 
SPECTRUM  OF  ARGON  IN 
THE  FAR  RED  REGION 
(ELECTRON  BEAM  PARALLEL 
TO  SPECTROGRAPH  SLIT) 
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FIG.  5.  9  A  SPECTRUM  OF  AN  ARC  HEATED  ARGON  PLASMA  JET 
-  Spectrogram  by  D.  Marsden,  UTIAS  - 
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FIG.  5.12  EXTENSION  OF  CALIBRATION  CURVES  TO  HIGHER  DENSITIES 


FIG.  5.12  EXTENSION  OF  CALIBRATION  CURVES  TO  HIGHER  DENSITIES 


f.  :  1.0  U  *  10 

p  s  5.44  torr  P„  *  9.0  torr 

O 


p*  ?  1 7  micron  s  4  micron 

Reo=2460  Re0  = 


o_ 

CN 


f0  *  0.425 
P0  =  2.79  torr 
Pa,  =  10  micron 

Re0=  840 


fo  =  0.425 
po  =  44.5  torr 
Poo  =  17.5  micron 
Re0=  4470 


_0> 

o 

u 

CO 


fo  =  0.12 

Pe  =  2.56  torr 
Poo  =  1  7  micron 
Re0=  533 


fo  =  0.12 

PQ  =29.7  torr 
P^  s  18.5 micron 
Re0=  2060 


FIG.  5.13 


FLOW  VISUALIZATION  OF  SOME  OF  THE  FREE  JETS 
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FIG.  5.14 


FLOW  VISUALIZATION  PHOTOGRAPH  OF  A  MACH  4.5 
ARGON  SHOCK  (FLOW  IS  FROM  LEFT  TO  RIGHT) 


FIG.  5.15 


DENSITY  PROFILE  THROUGH  A  MACH  4.5  ARGON  SHOCK 
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FIG.  5. 17 A  EXPERIMENTAL  ARRANGEMENT  FOR  THE  MASS-SPECTROMETER  TESTS 
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FIG,  5.17B 


SCHEMATIC  OF  "VEE  TUBE"  MASS-SPECTROMETER 
OF  THE  VEECO  MS-9AB  LEAK  DETECTOR 
(Reproduced  from  VEECO  Leak  Detector  Manual) 


FIG.  5.18  LEAK  DETECTOR  CALIBRATION  CURVE 
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FIG.  7.3 


PARTIAL  DENSITY  PROFILES  ALONG  THE  AXIS  OF  A  FREE 
JET  FOR  ReQ  =  227  AND  fQ  =  0.  21 
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FIG.  7.5 


PARTIAL  DENSITY  PROFILES  ALONG  THE  AXIS  OF  A  FREE 
JET  FOR  R«0  =  533  AND  f0  =  0.  12 
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FIG.  7.7 


PARTIAL  DENSITY  PROFILES  ALONG  THE  AXIS  OF 
A  FREE  JET  FOR  Re0  =  2060  AND  =  0.12 


FIG.  7.  9  MOLE  FRACTION  OF  ARGON  ALONG  THE  AXES  OF  VARIOUS  FREE  JETS 
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FIG.  7.  12  RADIAL  PARTIAL  DENSITY  PROFILES  FOR  A  FREE 
JET  OF  Ren  =  533  AND  fn  =  0.  12 


FIG.  7.13B  MOLE  FRACTION  OF  ARGON  THROUGHOUT  THE  FLOW  FIELD  OF 

A  FREE  JET  FOR  ReQ  =  533  AND  f  =0.12  (VIEW  IN  UPSTREAM 
DIRECTION) 


FIG.  7.14  RADIAL  PARTIAL  DENSITY  PROFILES  FOR  A  FREE  JET  OF 
Re  =  224  AND  f  =  0.  21 
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FLOW  VISUALIZATION  OF  A  SHOCK  WAVE  IN 
FRONT  OF  THE  SHOCK  HOLDER  (fQ  *  0. 49,  Mx  *  9) 
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FIG.  8.  3  PARTIAL  DENSITY  PROFILES  THROUGH  A  MACH  8.  7 
SHOCK  WAVE  (fQ  =  0.  11). 


FIG.  8.4 


PARTIAL  DENSITY  PROFILES  THROUGH  A  MACH  8.5 
SHOCK  WAVE  (fQ  =  0.  01) 


APPENDIX  A 


Analysis  of  the  Bottled  Gas  Mixtures 


The  majority  of  the  free  jet  and  shock  wave  experiments  we re 
performed  using  bottled  gas  mixtures.  After  the  first  free  jet  results  were 
analyzed,  it  became  clear  that  the  percentage  compositions  in  the  supply  cylin¬ 
ders  were  different  from  the  stated  values.  Consequently,  a  quantitative  analy 
sis  of  the  gas  mixtures  was  performed  by  making  use  of  the  experimental  setup 
available . 


Several  gas  flow  meters  (floating  ball  in  a  conical  tube) 
were  calibrated  for  both  argon  and  helium  at  2600  torr  using  as  a  standard 
an  "Elster,  Wet-Test  Model,  Research  Gas  Meter".  With  the  help  of  these 
calibrated  flow  raters  pure  helium  and  pure  argon  could  be  metered  into  the 
tunnel  from  separate  supplies,  so  that  any  percentage  composition  could 
easily  be  achieved  in  the  stagnation  chamber.  The  flow  through  the  flow 
meters  was  controlled  by  downstream  needle  valves,  across  which  most  of  the 
pressure  drop  occurred.  From  there  the  gases  flowed  into  a  low  pressure 
mixing  chamber,  which  was  connected  to  the  stagnation  chamber  by  a  tube  about 
4  m  long.  In  this  way  the  gases  were  mixed  at  low  pressure  (15  torr.), 
and  a  good  homogeneity  could  be  expected  before  the  gases  reached  the  nozzle. 

For  the  actual  calibration  run  a  5  mm  orifice  (nozzle)  was 
used,  the  stagnation  chamber  pressure  was  kept  constant  at  15  torr,  and  the 
100  m  amp  electron  beam  was  kept  at  a  fixed  position  (x  =  38  min)  on  the 
axis  of  the  frbe  jet  emerging  from  the  orifice.  The  helium  and  argon  signals 
from  the  optical  analyzer  were  recorded.  By  adjusting  the  amounts  of  helium 
and  argon  metered  in,  calibration  curves  of  these  photomuliplier  signals 
versus  mole  fraction  of  argon  in  the  source  were  obtained.  Because  of  the 
relatively  large  Reynolds  number  of  the  jet  (between  800  and  2000,  depending 
on  the  composition)  the  local  mole  fractions  of  argon,  calculated  from  the 
photomultiplier  signals,  differed  only  slightly  from  those  of  the  input  gas. 

These  calibration  curves  were  then  used  to  determine  the  mole 
fraction  of  argon  in  the  bottled  mixtures.  The  gas  cylinders  in  question 
were  alternately  connected  to  the  tube  leading  to  the  stagnation  chamber 
(via  a  pressure  regulator,  a  separate  flow  meter,  and  a  control  valve), 
and  electron  beam  intensity  measurements  were  made  under  exactly  the  same 
running  conditions  as  during  calibration.  The  helium  and  argon  signals 
obtained  from  the  optical  analyzer  then  determined  the  source  composition 
on  the  calibration  curves  (see  Fig.  A.l). 

The  nominal  values  of  the  gas  mixtures  and  the  calibrated 
values  were  as  follows: 


*0 

labelled 

calibrated 

cylinder  1 

0,05 

0.12 

cylinder  2 

0.25 

0.21 

cylinder  3 

0.50 

0.425 

cylinder  4 

0.05 

0.11 

cylinder  5 

0.50 

0.49 

Since  these  cylinders  had  been  standing  in  the  laboratory  for 
several  months  before  they  were  used,  it  is  highly  unlikely  that  the  gases  were 
not  thoroughly  mixed.  The  error  obviously  occurred  during  the  filling  opera¬ 
tion  of  the  cylinders  (the  composition  had  not  been  guaranteed  by  the  supplier) 


FIG.  A.  1  DETERMINATION  OF  THE  COMPOSITION  OF  THE  GAS  MIXTURES 
IN  THE  GAS  CYLINDERS 


«5 - - -  -  - 


APPENDIX  B 


Equilibrium  Value  of  fQ  in  the  Stagnation  Chamber 


The  equilibrium  value  of  fQ  in  the  stagnation  chamber  is  equal 
to  the  mole  fraction  of  argon  in  the  supply  cylinder  (or  to  the  volume  flow 
fraction  of  argon  metered  into  the  stagnation  chamber)  only,  if  the  integral 
of  0-fQ,  integrated  across  the  total  flow,  is  zero,  i.e.  if 


JJ  (0  -  fQ)  dS  =  0  ^  (B.l) 

S 

where  S  is  any  continuous  surface,  which  is  everywhere  normal  to  the  stream¬ 
lines,  and  which  crosses  all  streamlines  going  through  the  orifice.  Because 
of  continuity  of  the  mass  flow  of  each  species,  the  surface,  S,  need  not  be 
identical  with  the  sonic  surface  in  the  orifice. 


velocity  than  the  argon,  one  might  expect  the  above  integral  to  come  out 
negative;  the  positive  contribution  of  0  -  f  on  the  axis  being  outweighed 
by  the  negative  contributions  from  points  off  the  axis.  The  following 
argument,  however,  will  show  that  the  above  integral  vanishes  when  the  flow 
through  the  orifice  is  in  the  continuum  regime  (Re0  >  100). 


Consider  following  the  streamlines  backwards  through  the 
orifice  and  into  the  stagnation  chamber.  Many  orifice  diameters  (say,  more 
than  10)  upstream  of  the  orifice  the  streamlines  become  almost  straight,  and 
the  flow  becomes  spherically  symmetrical  and  closely  resembles  a  three-dimen¬ 
sional  sink  flow.  Now  let  the  surface,  S,  over  which  the  integral  is  to  be 
taken,  be  a  hemisphere*  of  radius,  R  (R  >  10D) ,  centered  around  the  orifice 
and  projecting  into  the  stagnation  chamber.  Then,  by  increasing  R,  the 
surface  area,  S,  increases  as  R^.  If  it  can,  therefore,  be  shown  that  at  every 
point  on  this  surface  the  quantity  0  -  f0  decreases  faster  t'^an  R"^,  as  the 
hemisphere  is  expanded,  then  the  integral,  //  (0  -  fQ)  dS,  mpst  vanish  as  R 
becomes  infinitely  large.  Since  the  surface,  S,  is  quite  arbitrary,  as  long 
as  it  is  normal  to  the  streamlines  and  extends  across  the  total  flow,  the 
surface  integral  must  also  vanish  for  any  permissible  finite  surface,  in 
particular  for  the  sonic  surface  in  the  orifice. 


*  or  a  surface  as  close  to  a  hemisphere  as  permitted  by  the  restrictions  on  S. 
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That  0  -  f  decreases  at  a  faster  rate  than  R"2  for  increasing 
R  can  be  shown  for  the  flow  close  to  the  axis  by  considering  Eq.  (l6)  of  Ref . 
8.  This  equation  is, 
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For  R  >  10D  the  flow  is  in  the  low  speed  subsonic  regime,  and 
the  gas  behaves  as  an  incompressible  fluid.  Therefore,  the  Mach  number  be¬ 
comes  simply  proportional  to  the  flow  velocity,  which,  according  to  the  con¬ 
tinuity  equation,  is  inversely  proportional  to  the  flow  area.  Hence,  as 
R  (or  x  -  00),  the  Mach  number  functions  occurring  in  Eq.  (B.2)  approach 

zero  in  the  following  way: 
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k 


M  = 


where  k  is  a  constant. 


R£ 


dM 

dM 

dx 

dR 

( dMN 

f 

4k2 

Vdx, 

>  = 

7 

1  ( 

'm) 

2 

1  V  — 

M  \ 

|  U.A  “ 

R 


4k 

F 


dR  =  -  ~ 

3  r3 


—00 


00 


Equation  (B.2)  then  becomes  of  the  following  form  for  R  -5 «  . 
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The  two  terms  in  Eq.  (B.2)  individually  decrease  as  R“3  as 
R  ->oo  .  Furthermore,  the  two  terms  become  exactly  equal  and  opposite  in 
sign  in  the  limit.  It  can,  therefore,  be  said,  that  f  -  0O  near  the  axis 
goes  to  zero  faster  than  R“3  for  R  becoming  large.  Because  of  the  spherical 
symmetry  of  the  flow  at  large  R,  it  is  unlikely  that  the  behaviour  of  0  -  fQ 
can  be  much  different  from  the  axial  behaviour  in  regions  away  from  the  axis. 


From  the  above  arguments  it  can  be  concluded  that  Eq.  (B.l)  is 
satisfied,  provided  the  Reynolds  number,  Re0,  is  large  (say,  Re0  >  100).  The 
latter  is  the  basic  assumption  underlying  Eq.  (B.2) 

This  means  that,  under  these  conditions,  the  number  flow  rates 
of  argon  atoms  and  helium  atoms  through  the  orifice  are  of  the  same  ratio  as 
the  corresponding  mole  fractions  in  the  stagnation  chamber.  Hence,  the  equili¬ 
brium  value  of  fQ  in  the  stagnation  chamber  is  equal  to  the  volumetric  argon 
flow  fraction  metered  into  the  tunnel,  which  in  turn  is  equal  to  the  mole  frac¬ 
tion  of  argon  in  the  supply  cylinder. 
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APPENDIX  C 

Relative  Thickness  of  the  Mach  Disc  in  a  Pare  Gas 


In  order  to  estimate  the  Mach  disc  thickness  in  terms  of  the 
jet  Reynolds  number  and  the  jet  pressure  ratio  across  the  orifice,  it  is  first 
necessary  to  find  the  state  of  the  gas  in  front  of  the  Mach  disc  as  a  function 
of  the  stagnation  chamber  conditions  and  the  ambient  pressure,  p  .  A  good 
estimate  of  the  Mach  number,  M-^,  in  front  of  the  Mach  disc  can  be  obtained 
by  making  the  following  two  assumptions: 

(l)  The  gas  flowing  along  the  axis  of  the  free  jet  expands  isentropically 
up  to  the  Mach  disc,  then  compresses  to  a  pressure  pP  >  p  through  a  normal 
shock  wave .  00 


(2)  The  pressure  behind  the  Mach  disc  is  simply  related  to  the  ambient 
pressure,  such  that  p0  =  €(?)p  ,  where  e{y)  is  independent  of  p  /p  . 

d.  00  o  00 

If  M^  were  known,  the  static  pressure  in  front  of  the  Mach 
disc  would  then  be  given  by  the  isentropic  relation. 
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The  pressure  p^  is  related  to  pg  by  the  Rankine-Hugoniot 
relations  for  a  normal  shock. 


p2  27M12  -7  +  1  €(7)  Poo 

_  _  7+I  =  p^“ 


Dividing  Eq.  (C.l)  by  Eq.  (C.2)  an  explicit  relation  is  obtained 
for  the  orifice  pressure  ratio  required  to  achieve  a  Mach  number,  M1?  in  front 
of  the  Mach  disc. 

Ill  p  7 A7-1) 

Po  e(7)  (1+2  Mi  )  (7+1)  /c  3s 

Poo  27M2  -7+1 

Alternately,  knowing  the  jet  pressure  ratio,  M^  can  be  obtain¬ 
ed  from  this  equation.  But  an  accurate  value  of  Mj_  can  also  be  found  inde¬ 
pendently  from  Wolff’s  calculations  (Ref.  13)  combined  with  the  semi- empirical 
Mach  disc  location,  given  by  Eq.  (l)  of  the  main  text.  The  M-j_  values  obtain¬ 
ed  from  Eq.  (C.3),  in  the  range  M]_  >  4,  can  be  brought  into  good  agreement 

with  those  given  by  the  characteristics  solution  of  Wolff  and  Eq.  (l),  by 
simply  assuming  that 

€(?)=—  .  (C.4) 

There  exists  no  direct  justification  for  this  assumption,  except  that  it 
facilitates  fitting  Eq.  (C.3)  to  the  calculations  of  Wolff. 

It  is  now  possible  to  combine  Eqs.  (C.3)  and  (-C.4)  with  Eq. 

(l)  in  order  to  obtain  a  semi-empirical  relation  for  the  axial  variation  of 
Mach  number  in  the  free  jet.  The  result  is, 
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which  agrees  quite  well  with  Wolff's  results  for  x/D  >  2. 


(C.5) 


Using  the  isentropic  relations  to  eliminate  M  from  Eq.  (C.5), 
the  state  of  the  gas  is  found,  in  terms  of  the  stagnation  state  for  any  point 
on  the  jet  axis.  For  x/D  >  2  one  may  thus  write 


(C.6) 


_  (y+  i)2  (y-  i) 

187  (x/D) 2 

Far  downstream  of  the  orifice  the  static  temperature  becomes 
very  much  smaller  than  T0,  and  Eq.  (C.6)  can  be  written  as 

_p_  _  /A1/7  =  ( l-Y^7”1)  -  (y+  i)2(y-i)  /c  , 

Po  \  V  187  (x/D)^ 

which  is  valid  only  at  large  x/D. 


Equation  (C.7)  predicts  that  at  high  Mach  numbers, many  ori¬ 
fice  diameters  away  from  the  orifice,  the  gas  density  decreases  inversely  as 
the  square  of  the  distance  from  the  orifice.  This  is  in  agreement  with 
Sherman's  source  flow  model  (Refs.  13  and  15),  which  predicts 


x0  =  location  of  hypothetical  source 
near  the  orifice, 


which  for  large  x/D  goes  to 


where  the  constant,  B,  is  chosen  such  as  to  fit  Wolff's  data.  As  will  be 
seen  from  the  following  table,  the  values  of  B  are  reasonably  well  predicted 
by  Eq.  (C.7) 
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(y+i)2(y-i)/i87 

5/3 

0.l6l 

0.158 

7/5 

0.0893 

0.0915 

9/7 

0.0615 

0.0645 
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The  Mach  disc  thickness  can  now  be  estimated  using  Eq.  (C.7)> 
if  it  can  be  assumed  that  xj^/D  is  large.  In  this  case  is  of  the  order 
of  10  or  larger.  Experimental  shock  thicknesses  in  this  range  are  about 
equal  to  4.5  mean  free  paths  ahead  of  the  shock  (see  e.g.  Ref.  34).  There¬ 
fore,  let  the  Mach  disc  thickness  be 


But, 


Therefore, 
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Using  now  Eq.  (C.7)  and  =  (T^T  )V  >  this  equation  can  be  re-written  as 

JL  =  1296  y3//2  R  -1  (  £l 

'yfzir  (7+l)2(r-l)  0  \  T0 

Substituting  Eq.  (l)  and  assuming  v  =  1/2  the  relative  Mach  disc  thickness 
becomes 
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